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Abstract 

A hardware model of the Gas Turbine Generator (GTG) in use on the US Navy's DDG- 
51 Class Destroyer is constructed for use as a lab apparatus at the Massachusetts Institute 
of Technology's Laboratory for Electromagnetic and Electronic Systems (LEES). The 
components of the hardware model include a 5 kilowatt three-phase generator; DC power 
supplies and motors that function as the prime mover; computer software to implement 
speed and voltage control; and an input-output interface board that passes measurement 
and controller signals to and from the software environment. 

A numerical Simulink model of the GTG is developed that provides speed response to a 
change in electrical loading. The GTG model takes into account basic physical 
characteristics of gas turbine generators and is tuned to provide a response similar to that 
of the destroyer's Allison 501-K34 GTG. An empirical open-loop model of the tabletop 
generator is also developed in Simulink and subsequently provided with closed-loop 
feedback control. Controller gains are adjusted such that the tabletop's Simulink model 
provides a response likened to the GTG model. 

Proportional and Integral (PI) control of the tabletop generator is implemented in the 
software environment. The tabletop generator's response to a certain electrical transient 
is compared to the GTG response predicted by the Simulink model. Recommendations to 
improve the response of the tabletop generator are made based on analysis of actual speed 
sensor noise. 
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Chapter 1 


Introduction 

Research being conducted at the Massachusetts Institute of Technology's Laboratory for 
Electromagnetic and Electronic Systems (LEES) includes shipboard applications of Non- 
Intrusive Load Monitoring (NILM) to reliably monitor and track diagnostic conditions of 
critical systems. The NILM's potential to disaggregate individual loads from bus current 
and voltage may provide future electrical protection systems with valuable information to 
protect critical loads [1], 

The Multi-Function Monitor (MFM) is currently employed to protect the Zonal Electrical 
Distribution (ZED) of modern ships. The MFM uses current sensors located on the main 
busses of the ship. This is considered a natural entry point for NILM to monitor multiple 
loads from a single point [1]. In addition to its current use on warships, ZED is a central 
concept in the architecture of the Next Generation Integrated Power System (NGIPS) [2]. 

Reliable technology is required to implement Power Distribution Modules (PDM) and 
Power Control (PCON) as called out in the NGIPS roadmap [2], The power monitoring 
being demonstrated at MIT is considered critical enabling technology for future PDM and 
PCON devices [1], 

This thesis describes the construction of a hardware model of the Gas Turbine Generator 
(GTG) in use on the US Navy's DDG-51 Class Destroyer. This hardware model of a 
ship's generator may be utilized in a scale model of a ship's electrical power distribution 
system. Subsequent simulations of shipboard electrical faults will allow for the 
advancement of power monitoring technologies. 

The equipment used in the hardware model is briefly described in the first two sections of 
Chapter 2. Section 2.3 details the characterization of the DC motor and evaluates the best 
suited gear ratio between the generator and prime mover. A potential gearbox is included 
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in the gear ratio discussion that may be utilized in lieu of existing timing belt pulleys. 
Chapter 3 provides an overview of gas turbines. This includes shipboard application of 
gas turbine generators; the basic design and classification of gas turbines; and discussion 
of gas turbine control and modeling. 

Simulink models of the GTG and the tabletop generator are developed in Chapter 4. The 
GTG model follows from the control and modeling discussed in Section 3.3. The GTG 
model is initially tuned to meet response specifications. Minor adjustments are then 
made to key parameters in the model to match actual response data. 

Section 4.2 discusses the empirical development of the open-loop model for the tabletop 
generator. This system identification is performed in two distinct steps. First, a Simulink 
model simulates the response to a change in electrical loading while terminal voltage to 
the DC motor is held constant. This captures the self regulating aspect of a DC motor as 
discussed in Section 2.3. Building on this model, the second step incorporates the effect 
that a change in terminal voltage has on speed. In this open-loop test, a change in 
terminal voltage drives the speed response with the generator's electric load held 
constant. 

Following the open-loop modeling of Section 4.2, closed-loop feedback control of the 
tabletop generator is modeled in Section 4.3. The closed-loop model is tuned to provide 
a response similar to the GTG model for a specific change in electrical loading. The 
closed-loop model is then evaluated for the following: sensitivity to a key parameter that 
showed deviation in open-loop modeling; exclusion of estimated power losses; and 
variation in step changes of the electric load. 

A physical description of the hardware model is provided in Chapter 5 along with 
response testing results. Recommendations to improve the response of the tabletop 
generator are made based on analysis of actual speed sensor noise. Conclusions and 
future work are discussed in Chapter 6. 
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Implementation of the control scheme in the software environment was conducted by 
Vanessa Esch while working on an Undergraduate Advanced Project. Signal 
conditioning of speed and voltage measurements was developed by Jacob Osterberg 
while working through the Undergraduate Research Opportunity Program. Their work is 
shown in its entirety in the Appendices. Though some information is duplicated and 
conflicts may exist, no attempt was made to modify or summarize their reports. This 
preserves the insight of three observers involved on one project. 
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Chapter 2 


Initial Design Considerations 

The equipment used in the hardware model is briefly described in the first two sections of 
Chapter 2. Section 2.3 details the characterization of the DC motor and evaluates the best 
suited gear ratio between the generator and prime mover. A potential gearbox is included 
in the gear ratio discussion that may be utilized in lieu of existing timing belt pulleys. 

2.1 Generator 

Obtaining a reasonably sized three phase AC generator was the starting point for creating 
the tabletop model of a shipboard generator. Some initial effort was put forth in utilizing 
a large truck alternator which produces a three phase current prior to its rectification for 
use in the truck's DC power system. However, given that alternators are intended to 
produce DC power and are thus not restricted in their speed; operation of the alternator at 
a speed to provide a frequency of 60 Hz results in an undesirably low power output. 
Alternator conversion kits exist that provide reasonable 60 Hz power output. Conversion 
of an alternator requires additional windings on the armature to provide more power at a 
slower speed. This was deemed impractical, removing the alternator from further 
considerations. 

Three phase generators are typically associated with large industrial applications and are 
not readily available in smaller power ratings. The search for a smaller three phase 
generator resulted in the purchase of the 5 KW machine shown in Figure 1. Generator 
nameplate data is shown in Figure 2. These generators are primarily intended for use in 
remote or otherwise underserved areas where electrical power is either not available or 
unreliable. Typically, these machines would be driven by a tractor's Power Take-Off 
(PTO) drive or any other available combustion engine and utilized where exacting control 
of frequency is not required. Additional specification is included in Appendix C and 
purchase information is listed in Appendix D. 
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Figure 1: Generator Head 
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Figure 2: Generator Nameplate 
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2.2 Prime Mover 


A direct current motor was chosen to act as the generator's prime mover due to its 
inherent controllability as discussed in Section 2.3. Two permanent magnet DC motors 
chosen from lab inventory are shown coupled together in Figure 3. Each motor is rated at 
2 HP as indicated on the nameplate shown in Figure 4. The motors have a drive shaft 
that extends outward from either side of the motor housing which allows the motor to be 
centered between two driven loads such as on a treadmill. Rigidly coupling two motors 
together and connecting them in series electrically provides a rated capacity of 4 HP. 
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Figure 4: Motor Nameplate 


Remote output control of a DC power supply was required to implement speed control of 
the prime mover. The remote output control feature of a power supply allows a small 
voltage input, such as from a PC, to direct the voltage output from the power supply to 
the motor. Two XHR 1000 Watt Series programmable DC power supplies were chosen 
from lab inventory to meet this requirement. The power supplies were placed in a 
parallel master-slave configuration to provide a rated output of 150 volts and 14 amps. 


2.3 DC Motor Characteristics 


The circuit diagram for a permanent magnet DC motor is shown in Figure 5. The 
inductance of the armature windings is typically neglected in the study of simple DC 
motors. Kirchhoff s Voltage Law leads directly to the basic equation of a DC motor as 
shown in Equation 1 [3]. This is the starting point to develop the equations necessary to 
characterize the DC motor. The equations that follow are utilized in the identification of 
the best suited gear ratio at the end of this section. The DC motor characteristics and 
equations are also used in Sections 4.2 and 4.3 where Simulink modeling of the tabletop 
generator includes the mechanical loading of the motor. 
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Figure 5: DC Motor Circuit Diagram 

Voltage applied to the terminals of the DC motor produces a current on the armature 
windings. This current is within the influence of a magnetic field produced from either 
permanently installed magnets or field windings. The current carrying conductor within 
a magnetic field satisfies the requirements for motor action, producing a force that acts 
orthogonal to the current flow. A commutator maintains the orientation of the armature 
current to ensure this force results in motor rotation. Upon rotation of the motor, the 
conductor travelling within a magnetic field satisfies the requirements for generator 
action, producing a voltage that counters the applied voltage [3]. 

v Term =i M RM+v c (i) 

Where: V Ter m = terminal voltage 

Im = armature current 
R m = armature resistance 
Vc = counter voltage (back electromotive force) 

If the magnetic field is maintained constant through the use of permanent magnets or by a 
constant current applied to field windings, then the counter voltage of a DC motor is a 
linear function of its speed as shown in Equation 2 [3]. 
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V c = Kco (2) 


Where: K = motor constant 

CO = rotational speed (radians/sec) 

The value of the motor constant was determined by mechanically joining two motors as 
shown previously and applying a range of voltages to the terminals of one of the two 
motors. With one motor driving the other, the counter voltage of the second at various 
speeds may be measured directly at its terminals with no current on the armature. Linear 
regression of the counter voltage data plotted as a function of rotational speed results in a 
straight line with a slope equal to the motor constant. The result for one of the motor 
constants is shown in Figure 6 to be 0.215 volt*sec. The second motor tested produced a 
similar value with a motor constant of 0.217 volt*sec. 


Motor Constant 



Figure 6: Motor Constant 
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Power in a simple DC circuit is the product of current and voltage. The voltage drop 
associated with the armature resistance shown in Equation 1 results in some power lost to 
heat in the motor. Neglecting frictional heating, the useful mechanical power of the 
motor (Pm) becomes the product of armature current and counter voltage and may be 
written as: 


P m =KI m co (3) 

This result may also be obtained by considering the armature torque that is produced 
from the interaction of the field and armature flux [3]. This torque is a result of the motor 
action force mentioned earlier acting at some distance from the center of rotation. Given 
the constant field flux of the permanent magnet motor, the torque of the motor (T M ) may 
be written as: 


T m =KI m (4) 

Mechanical power in its simplest form arises from force acting over some distance for 
some period of time. For rotating machinery, power is the product of torque and 
rotational speed as shown in Equation 5. 

Pm=T m cq=KI m O) (5) 

Direct current motors have long been utilized in applications requiring speed control. 

The linear nature of the equations presented above attest to the controllability and 
inherent stability of the DC motor [3]. For example, increasing the mechanical loading 
on the motor will initially slow the machine and decrease its counter voltage. Since 
armature resistance is typically very small, any decrease in counter voltage results in a 
pronounced increase in armature current. This balances out until a new steady state is 
reached where the power of the motor equals the mechanical loading at some lowered 
speed with an increased armature current. In order for the DC motor to return to its initial 
speed, some operator or automatic control function must occur to raise terminal voltage. 
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After obtaining the motor constant, a torque-speed curve was developed to further 
classify the DC motor and provide insight into acceptable motor speeds. Combining 
Equations 1 and 2, solving for current, and substitution into Equation 4 results in the 
following: 


T„=-^V T , rm -^-o> (6) 

'\V7 

Motor resistance was determined at the same time as the motor constant with one motor 
driving the other. For motor resistance, the counter voltage developed on the slowly 
driven motor was measured at a specific speed. Shorting the driven motor's terminals 
through a DC meter and readjusting speed provided the current flow corresponding to the 
previously measured voltage. Dividing voltage by current resulted in a motor resistance 
of 0.957 and 0.912 ohms for the two motors tested. With motor constant (K) and motor 
resistance (Rm) known, Equation 6 was utilized to show how motor torque changes with 
speed for a given applied terminal voltage. The resulting torque-speed curve for both 
motors mechanically joined and operating in series electrically is shown in Figure 7. 
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Torque-Speed for 2 Motors in Series 



Figure 7: Motor Torque-Speed 

The two motors operating in series have a rating of 4800 rpm at 240 volts; however, the 
power supply is limited to 150 volts. The torque-speed curve was plotted at 100, 125, 
and 150 volts to allow for evaluation of motor speed given this constraint. From 
Equations 1 and 2, the motor speed corresponding to a power supply of 150 volts and 14 
amps was determined to be 2737 rpm. 


co = 


V T . 


I M^M 
K 


(7) 


Thus, the best suited motor speed would be 2737 rpm to properly deliver power to the 
prime mover from the power supply. The torque-speed curve indicates that the motor 
will operate satisfactorily at this speed between 125 and 150 volts. Interpolation of 
unloaded speed data collected at the same time of motor constant testing indicates the 
unloaded motors should operate very near 125 volts at 2737 lpm drawing 0.7 amps. 
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Given the generator must rotate at 1800 rpm to produce power at 60 Hz; the best suited 
gear ratio is 1.52. Any deviation from this gear ratio will limit the maximum output of 
the power supply. From Equation 7, in order for the motor to operate at a higher speed, 
given that the terminal voltage is limited in value, the motor current must decrease. Thus, 
for a larger gear ratio, the power supply will become voltage limited, reaching 150 volts 
at a current less than the rated peak. Conversely, a smaller gear ratio will result in a 
current limited power supply where 14 amps is reached prior to 150 volts. 

For the purposes of this work, essentially a thorough feasibility study, the motor and 
generator were mounted aside one another on a rigid platform. The motor was 
mechanically coupled to the generator via L Series timing belt pulleys and a 1/2 inch 
timing belt. This is later shown in Chapter 5 with the physical description of the 
hardware model. Given the initial purchase of a generator pulley with 32 teeth from 
preliminary work involving one 1.5 HP motor, the readily available gear ratios were 
either 1.6 or 1.4545 corresponding to a new motor pulley with 20 or 22 teeth respectfully. 

By fixing either terminal voltage at 150 volts for ratios larger than 1.52 or current at 14 
amps for smaller ratios, Equation 7 can be used to show how various gear ratios limit the 
use of the power supply. Table 1 compares the two ratios above for this initial design and 
includes two variations. The first variation is the use of HTD Series pulleys available 
with 34 and 22 teeth corresponding to a gear ratio of 1.5454. The second variation is an 
angled Howse Gearbox with a ratio of 1.47 that is included for future consideration. 

The Howse Gearbox is designed for a tractor's PTO driven M60 Rough Cut Mower. A 
tractor's PTO drive operates at 540 rpm to drive a wide assortment of attachments. As 
such, the intended speed range of this gearbox is approximately 1/5 that of the tabletop 
motor driven generator. However, it is rated at 55 HP, more than 10 times the power 
driving the generator, and is designed to withstand significant backlash associated with 
objects hit by a large tractor pulled mower. Given the gearbox's splined shafting and its 
angled design, it is expected that a fair amount of machining and alignment work would 
be required for its incorporation. 
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Specifications and ordering information is included in Appendix D for both the hardware 
utilized to achieve a gear ratio of 1.4545 and the Howse Gearbox that may be 
incorporated in future work. 



Gear Ratio 

Motor Speed 
(rpm) 

Max Power Supply 
Output (watts) 

L Series (32/20) 

1.6 

2880 

1581 

HTD Series (34/22) 

1.5454 

2782 

1938 

Best Suited Ratio and Speed 

1.52 

2736 

2100 

Howse Gearbox 

1.47 

2646 

2042 

L Series (32/22) 

1.4545 

2618 

2025 


Table 1: Gear Ratio Comparison 


The gear ratios shown in Table 1 that are less than the best suited ratio of 1.52 represent 
the current limited power supply as discussed previously. This table identifies that it is 
more advantageous to operate below the best suited speed with regard to utilization of the 
given power supply. Operation at a slower speed is also beneficial from a torque-speed 
standpoint as shown in Figure 7. As such, the L Series 22 tooth timing belt pulley was 
installed on the motor shaft to drive the 32 tooth pulley on the generator. Table 1 also 
shows that the Howse Gearbox has a well suited gear ratio for this application. 
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Chapter 3 


Gas Turbine Generators 

Chapter 3 provides an overview of gas turbines. This includes shipboard application of 
gas turbine generators; the basic design and classification of gas turbines; and discussion 
of gas turbine control and modeling. The modeling of a GTG control system discussed in 
Section 3.3 leads directly to the Simulink model of the Allison 501-K34 GTG shown in 
Section 4.1. 

3.1 Shipboard Application 

As will be discussed throughout Sections 3.2 and 3.3, single-shaft gas turbines are 
inherently better suited for constant speed applications and are typically preferred for 
production of electrical power when compared to two-shaft gas turbines. The Allison 
501-K34 GTG installed on US Navy DDG-51 Class Destroyers is a single-shaft gas 
turbine with a compression ratio of 10.5 operating at a speed of 14,340 rpm. It is rated to 
provide a usable output of 4328 HP without bleed air in use and 3523 HP while supplying 
2.37 lb/sec of bleed air [4], Bleed air is air that may be removed from the compressor to 
provide compressed air to various auxiliary systems throughout the ship. 

Given the specifications above, the Allison 501-K34 would be classified as a light 
industrial gas turbine engine with a high power-to-weight ratio as described at the 
beginning of Section 3.2. The power ratings correspond to 3227 KW and 2627 KW 
respectfully. This variation in power rating, dependant on the utilization of bleed air, 
accounts for discrepancies seen in early literature review where some documents rate the 
Allison 501-K34 at 3 MW while others identify its rating as 2.5 MW. 

The governing system of the shipboard GTG is required to provide isochronous control 
while allowing for load sharing with other generators operating in parallel with similar 
governing systems [4]. This is the normal mode of operation with two or three generators 
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providing power to the ship. An additional droop control mode is included to allow the 
shipboard GTG to load share with pier side shore-power. 

In addition to providing isochronous governing with load sharing capability, the Allison 
501-K34 Specification requires the governing system to handle full power load 
excursions. That is, without bleed air in use, the GTG must withstand the application or 
removal of its full rated load without exceeding a 2% deviation from its base speed and 
must recover within 1.5 seconds to stay within 1% of its base speed [4], The strict speed 
control placed on the Allison 501-K34 driven generator along with its high power-to- 
weight ratio attest to the demanding requirements of naval warship machinery. 

3.2 Basic Design 

Gas turbines are commonly classified based on their intended use and typically fall into 
one of the following categories: aerospace, light industrial, or heavy industrial. Light 
industrial gas turbines have high power-to-weight ratios obtained through operation at 
high combustion and exhaust temperatures along with relatively high compression ratios 
that typically exceed 7. Structural weight is minimized which generally implies more 
frequent maintenance. To offset this, light industrial turbines are designed as a single 
module for easy removal and installation. Thus, a malfunctioning unit may be quickly 
swapped for a functional gas turbine and sent out for repair. Light industrial gas turbines 
are limited in their power output to about 10 MW [5], 

Gas turbines are also classified based on their physical construction. In a single-shaft gas 
turbine, the compressor and power turbine share a common shaft. Two-shaft gas turbines 
have a compressor turbine that drives the compressor located on a common shaft while 
the power turbine resides on a separate shaft to drive the load. The basic arrangement of 
a single and a two-shaft gas turbine are shown below [5]. 


23 




Figure 8: Basic Arrangement of a Single-Shaft Gas Turbine 



Figure 9: Basic Arrangement of a Two-Shaft Gas Turbine 


Two-shaft gas turbines allow both the compressor and the power turbine to operate at 
their most efficient speeds, allowing for better thermodynamic performance. However, 
two-shaft gas turbines are slower to respond and are not suitable to drive a generator that 
must operate in synchronism with other generators [5]. 
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The basic operation of the gas turbine is best described by the Ideal Brayton Cycle that 
consists of two isobaric and two isentropic processes. Figure 10 shows the applicable 
Pressure-Volume (P-V) and Temperature-Entropy (T-S) diagrams that describe the 
thermodynamics of a single-shaft gas turbine [61. 
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Figure 10: Ideal Brayton Cycle 


Referring to Figure 8 and Figure 10, the compressor works on the fluid from 1 to 2, 
increasing pressure and temperature while decreasing volume at constant entropy. 
Combustion occurs from 2 to 3, increasing temperature, entropy, and volume at a 
constant pressure. The compressed and heated fluid expands through the power turbine 
to perform work from 3 to 4 at constant entropy with an associated decrease in pressure 
and temperature. Heat is removed from the system from 4 to 1 via the exhaust. 
Regenerative cycles are often included to capture some of the heat lost through the 
exhaust [6]. 

3.3 Basic Control 

Control systems ensure a system produces a desired output and are comprised of control 
loops that may be categorized as either open or closed. In an open loop control system, 
the input is independent of the output and there usually exists an offset between the 
output and set point. In the closed loop system, the offset from the open loop response is 
input to the closed loop controller to moderate the output by adjusting the input. The 
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closed loop control systems predominately found in use incorporate some combination of 
proportional (P), integral (I), and derivative (D) action [7]. 

Proportional control produces a signal proportional to the error between the measured 
system output and the desired set point. Proportional control alone will result in a steady- 
state error known as proportional offset if the set point is adjusted. The offset may be 
alleviated either through operator action (manual reset) or by the inclusion of integral 
action in the closed loop controller [7], 

Integral action occurs as a result of the error being integrated continuously and 
effectively provides an automatic reset that removes proportional offset. One drawback 
to inclusion of integral control is the possibility of integral wind-up that may occur if 
conditions are met such that controller action no longer affects the system output. This 
may result from actuator saturation where some physical limitation is reached that 
prevents further response. The resulting steady error will be summed up by the integral 
action and force unnecessary control when the system re-enters a controllable range [7], 

Derivative control looks at the rate of change of the error and is useful to improve 
transient response when the system is inherently very slow. Derivative action will not 
counter proportional offset or integral wind-up given these occur with a non-changing 
error. Derivative control is often omitted in gas turbine control systems. A controller 
utilizing all three terms (PID) is represented below [7]. 

OP = Er-Kc + —\Erdt + Kc ■ Td (g) 

Ti J dt 


Where: OP = controller output 

Er = error between process output and set point 
Kc = controller gain 
Ti = integral time 
Td = derivative time constant 
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Gas turbine control systems must also protect the engine from exceeding design limits. 
This is accomplished with the incorporation of signal selection where additional error 
signals are generated from the comparison of operating limit set points with current 
values. Two examples of operating limits that preclude engine damage are Exhaust Gas 
Temperature (EGT) limit and power turbine speed limit. Signal selection may also be 
employed to avoid operation near critical speeds where a significant increase in vibration 
occurs, typically 50 to 70% of design speed for a gas turbine [7], 

Single-shaft gas turbines are most suited for fixed speed applications with a resistance to 
over-speeding resulting from the high power requirements of the compressor. They are 
typically equipped with Variable Inlet Guide Vanes (VIGV) on the compressor to reduce 
starting power requirements by reducing airflow through the compressor. The position of 
the variable guide vanes are adjusted based on EGT. Two schemes exist where one 
maintains EGT below the operational limit and the other maintains it at the limit during 
reduced power output. The second is typically employed with combined cycle plants 
where the hot exhaust gas from the gas turbine is utilized to operate or supplement a 
steam cycle power plant. In either case, another control scheme is required in addition to 
that maintaining speed or power [7]. 

The open-loop or uncontrolled speed-torque characteristic of a gas turbine is not suitable 
for power generation. Without closed-loop feedback control, the initial speed response to 
a change in electrical loading is predominately a function of the generator's inertia. A 
closed-loop gas turbine governor measures and amplifies deviation in speed, sending a 
response signal to the fuel valve that regulates gas turbine power. At this point, there is 
an inherent variation between the response of single-shaft and two-shaft gas turbines. 
Two-shaft gas turbines have an additional finite delay in that the compressor responds 
prior to the power turbine. As such, single-shaft gas turbines are considered to have 
superior speed performance, able to accept sudden changes in electrical power with less 
deviation and faster recovery [5]. 
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A generator's output frequency is determined by its speed. In all power systems, 
frequency deviation must be minimized for incremental changes in power demand, some 
of which could be quite large. There are two methodologies used predominately to 
control the speed of gas turbine generators. These are droop governing and isochronous 
governing [5]. 

Droop governing, also referred to as proportional control, allows a drop in speed to occur 
as a result of increased electrical loading. Droop governing is commonly used in power 
systems since it provides a simple and reasonably accurate load sharing capability 
amongst a group of generators [5]. This method reflects the proportional offset discussed 
earlier that occurs from proportional only control. 

Isochronous governing, also referred to as integral control, drives the steady state speed 
error to zero, resulting in a constant frequency. Accurate power sharing and constant 
speed control involves load measurement of each generator, power system frequency 
measurement, and a control sub-system that drives power mismatches to zero [5]. For a 
single generator, isochronous governing may be implemented through proportional and 
integral control with a relatively large gain placed on the integral action. Accurate load 
sharing amongst two or more isochronous generators is beyond the scope of this work. 

A basic block diagram is shown in Figure 11 that represents the elements of the equation 
of motion that are most associated with gas turbine control. Rotational friction and 
windage are often ignored since they have little influence on the performance of the 
control system. The complexity of a gas turbine control block diagram is dependent on 
available data and the nature of the study. Simplifications are usually acceptable, for 
example, engine protection sub-systems are not required for basic response analysis [5]. 
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Figure 11: Basic GTG Control Block Diagram 


A model representing the droop governing of a single-shaft gas turbine is shown in 
Figure 12. Electrical power is the actual load demand on the power turbine shaft and 
mechanical power is the useful power output of the gas turbine. Electrical power is 
subtracted from mechanical power such that an increase in electrical loading initially 
results in a decrease in speed as determined by the inertia of the system. The system 
inertia includes the gas turbine, couplings, gearbox, and generator. The term representing 
this inertia as shown in the model is defined as [51: 


Gh = 


1 

2 B 


(9) 


Where: Gh = inverted inertia term 

H = inertia constant (seconds) 

The inertia constant (H) is typically used in electrical engineering to convert the moment 
of inertia of the rotating system to a base of electrical volt-amps. It was developed 
specifically for solving differential equations describing generator shaft dynamics and is 
defined as the energy stored in the rotating mass divided by the volt-amp rating of the 
generator [5]. 
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Figure 12: Droop Governing of a Single-Shaft Gas Turbine 


All values in the model are normalized to a per unit system which allows for evaluation 
of machines with widely varying capacities. Power is normalized to the rating of the 
generator and speed to that which corresponds to a frequency output of 60 hertz. The 
speed change from the inertia transfer function is summed with the per unit base speed. 
The resulting turbine speed is then subtracted from the setpoint to provide an error signal 
to the droop governor. Governor lags 1 and 2 represent delays inherently present in 
electronic circuits and contain the derivative damping gain Kg2 which is often adjustable. 
A compensation circuit, used in some controllers to improve speed response, is 
represented by the governor lead-lag block [5]. 


The output signal from the controller directs the fuel valve position which adds an 
additional delay associated with the valve stem response. The fuel valve has physical 
limits represented in a saturation block that prevents the model from demanding a greater 
change in fuel than what is achievable. An artificial negative value for the minimum fuel 
flow represents the fact that the valve is initially open to about 15% of its travel while the 
gas turbine generator is operating at speed with no electrical loading [5]. 
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There is a finite burning rate the associated with complete combustion of the fuel as it 
travels through the combustion chamber. Additionally, there is a delay as the heated gas 
transfers energy through the power turbine to produce usable mechanical power. The 
turbine lead-lag approximates this conversion of a change in fuel flow to a change in 
power [5]. 


A range of typical values found in the gas, oil, and petrochemical industry is provided in 
Table 2. The values shown in this table are based on per unit modeling of the gas turbine 
control system [51. Applicable parameters are utilized in the modeling of the Allison 
501-K34 GTG as presented in Section 4.1. Values shown here serve as the starting point 
for the Simulink model of a shipboard GTG. 


Parameter 

Low 

Typical 

High 

Gh 

0.25 

0.33 

0.42 

Kdg 

0.02 

0.04 

0.08 

Kgl 

1.0 

1.0 

1.0 

Tgl 

0.015 

0.01 

0.05 

Kg2 

10.0 

20.0 

40.0 

Tg2 

0.02 

0.04 

0.15 

Kg3 

1.0 

1.0 

1.0 

Tg3 

0.25 

0.50 

0.75 

Tg4 

1.0 

1.50 

1.75 

Tf 1 

0.01 

0.02 

0.05 

fmax 

1.2 

1.35 

1.5 

fmin 

-0.2 

-0.15 

0.0 

Ttl 

0.3 

0.6 

0.9 

Tt2 

1.2 

1.4 

2.0 


Table 2: Typical GTG Parameters for the Gas, Oil, and Petrochemical Industry 
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Chapter 4 


Model Development 

Simulink models of the GTG and the tabletop generator are developed in Chapter 4. The 
GTG model follows from the control and modeling discussed in Section 3.3. The GTG 
model is initially tuned to meet response specifications. Minor adjustments are then 
made to key parameters in the model to match actual response data. 

Section 4.2 discusses the empirical development of the open-loop model for the tabletop 
generator. This system identification is performed in two distinct steps. First, a Simulink 
model simulates the response to a change in electrical loading while terminal voltage to 
the DC motor is held constant. This captures the self regulating aspect of a DC motor as 
discussed in Section 2.3. Building on this model, the second step incorporates the effect 
that a change in terminal voltage has on speed. In this open-loop test, a change in 
terminal voltage drives the speed response with the generator's electric load held 
constant. 

Following the open-loop modeling of Section 4.2, closed-loop feedback control of the 
tabletop generator is modeled in Section 4.3. The closed-loop model is tuned to provide 
a response similar to the GTG model for a specific change in electric loading. The 
closed-loop model is then evaluated for the following: sensitivity to a key parameter that 
showed deviation in open-loop modeling; exclusion of estimated power losses; and 
variation in electric loading. 

4.1 Numerical Simulink Model of Allison 501-K34 GTG 

The numerical Simulink model developed to simulate the response of the shipboard 
generator is shown in Figure 13. This model is similar to that shown in Figure 12 with 
the notable difference that the droop governing control blocks are replaced with a PI 
controller. As discussed in Section 3.3, a large gain placed on the integral action of a PI 
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controller invokes isochronous control. Thus, the modeled controller represents the 
constant speed governor of the GTG during normal shipboard operations. 



Figure 13: Simulink Model of Allison 501-K34 


As with the droop governed model, per unit values are utilized to allow for comparison 
with machines of significantly different ratings. Aside from the governor, the core 
characteristics and signal flows are the same as that discussed in Section 3.3. The inertia 
term (Gh), as defined in Equation 9, was estimated from the inertia of the gas turbine and 
generator. The mass moment of inertia for the gas turbine was taken from the Allison 
501-K34 Specification [4], The generator's inertia was obtained via correspondence with 
Naval Surface Warfare Center, Carderock Division (NSWCCD) [8]. 


Initially, typical values from Table 2 were used in the model to test its response to a step 
increase in electrical power from 0 to 1 per unit, representing the instantaneous 
application of the machine's full rated load from 0 to 3000 KW. Controller gains were 
increased and time constants adjusted until the speed response met the GTG Specification 
criteria discussed in Section 3.1. Table 3 shows typical ranges of the applicable 
parameters taken from Table 2 along with the values utilized in this model to meet speed 
response criteria. As with Table 2, the parameters are based on per unit modeling of the 
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speed control of a GTG typically found in the gas, oil, and petrochemical industry. They 
are included here to provide some validation of the shipboard GTG model. 


Parameter 

Typical Values 

Simulink Model 

Gh 

0.25 to 0.42 

0.55 

Kp 

- 

40 

Ki 

- 

40 

Tf 1 

0.01 to 0.05 

0.01 

fmax 

1.2 to 1.5 

1.35 

fmin 

-0.2 to 0 

-0.15 

Ttl 

0.3 to 0.9 

1.2 

Tt2 

1.2 to 2.0 

0.9 


Table 3: Initial GTG Simulink Model Parameters 


Although the governor of the shipboard GTG is not precisely modeled and it is not 
known whether the turbine lead and lag time constants are truly representative; the model 
does include the fundamental dynamic aspects of gas turbine generators and values 
utilized in the model are typical or near typical as shown in Table 3. 

A larger value of the inertia term (Gh) represents a machine that has less stored rotational 
energy given the same power rating. Given space and weight limitations of naval vessels, 
shipboard machinery often has larger than normal power densities and it is not surprising 
that this value is outside the typical range. However, the inertia of the gearbox was not 
specifically included with the value obtained for the generator and while the inertia of the 
generator (3550 lb ft') clearly dominates when compared to the gas turbine (61 lb ft'); the 
gearbox may represent a reasonable percentage of the system's inertia. If the inertia of 
the gearbox is not included in the generator's value, and assuming the gearbox accounts 
for 10% of the system's inertia, then the value for Gh becomes 0.50 which is closer to but 
still outside the typical range. 
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Given that exact model parameters are unknown and the intent was to provide a response 
similar to the GTG, the model was deemed acceptable given its ability to meet response 
criteria as shown in Figure 14. 


GTG Model Speed Response 



Figure 14: Initial GTG Model, Response to Rated Load Step Increase 


Following the development of the model above, actual full load power excursion data 
was received from NSWCCD [8]. As tested, the Allison 501-K34 response to removal 
and application of its full rated load is shown in Figure 15. The associated electrical 
transient is shown in Figure 16. 
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Figure 15: Actual GTG Response to Full Rated Power Excursions 



Figure 16: Actual Full Rated Electrical Power Transient 









Utilizing the provided data, the actual full rated electrical power transient was input to the 
GTG model vice an idealized step input. The comparison of the model's speed response 
with actual data is shown in Figure 17. 



Figure 17: Initial GTG Model Comparison 


As seen in Figure 17, the actual response exceeds the 2% deviation initially suggested by 
the 1987 Model Specification. This is most likely attributable to variations found on 
ships constructed later in the class (DDG-91 and beyond). Modifications were made to 
the model parameters to more closely approximate the actual response. The model is 
most sensitive to small changes in the turbine lead and lag time constants Ttl and Tt2 
respectfully. Some improvement was also noted utilizing a smaller value for the inertia 
term (Gh). The inertia term was reduced in value to that mentioned earlier, accounting 
for the inertia of the gearbox, with further reduction showing little affect. No changes 
were made to the modeled controller gains or fuel valve dynamics. The modified 
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parameters are shown in Table 4. The response of the finalized model of the GTG as 
compared to the actual response data is presented in Figure 18. Although the model does 
not correctly identify the overshoot seen in the actual data, it does adequately describe the 
peak response and settling time. This model becomes the basis for the tabletop 
generator's empirical model developed in Section 4.2. 


Parameter 

Typical Values 

Initial Value 

Final Value 

Gh 

0.25 to 0.42 

0.55 

0.50 

Kp 

- 

40 

40 

Ki 

- 

40 

40 

Tf 1 

0.01 to 0.05 

0.01 

0.01 

fmax 

1.2 to 1.5 

1.35 

1.35 

fmin 

-0.2 to 0 

-0.15 

-0.15 

Ttl 

0.3 to 0.9 

1.2 

0.96 

Tt2 

1.2 to 2.0 

0.9 

1.14 


Table 4: Modified GTG Simulink Model Parameters 
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Modified Model Comparison with Actual Data 



Figure 18: Finalized GTG Simulink Model 


4.2 Empirical Open-Loop Modeling of Tabletop Generator 

Open-Loop response testing was performed for the following conditions: increased 
electrical loading from 0.045 to 0.135 per unit, decreased electrical loading from 0.135 to 
0.045 per unit, and increased terminal voltage by 5.2 volts while supplying 0.045 per unit 
power near operational speed and output voltage. The per unit changes in power 
correspond to the readily available 225 to 675 watt step loading of the tabletop generator 
rated at 5 KW. The empirical model of the tabletop generator's speed response to a 
change in electrical load as developed in Simulink is shown in Figure 19. 
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Figure 19: Open-Loop Response Model to a Change in Loading 


The model was designed in a fashion similar to the GTG model developed in Section 4.1 
where a mismatch between electrical and mechanical power results in an initial change in 
speed as determined primarily by the inertia of the machine. The tabletop generator 
presents a fair amount of vibration throughout all speeds in addition to pronounced 
vibration that occurs at a number of critical speeds. The sideward loading associated 
with the prime mover pulley arrangement may exacerbate this condition. The power loss 
approximates this loss of energy and allows for a more direct observation of motor power 
within the model. The lagging transfer function describes how a change in speed affects 
the motor's power output. This self regulation, within the uncontrolled or open-loop 
model, attests to the DC motor's inherent stability as discussed in Section 2.3. 


Figure 19 shows initial values for power and speed immediately prior to a step increase in 
electrical loading. Initial motor power was determined from Equation 3, utilizing 
measured values for motor current and speed. Initial estimates for the prime mover's time 
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constant and gain (Tpm and Kpm) were made utilizing techniques most valid for a 
system dominated by a primary lag. The time constant was taken as the time required for 
the response to reach 63% of the total change in speed. The gain was estimated as shown 
below. Figure 20 provides a graphical representation [9]. 

Kpm = (10) 

A u 

Where: Aco = completed response change in speed 

A u = straight line change in speed 


Speed 


Time 

Figure 20: Initial Estimation of Prime Mover Gain 

Using these estimated values, the inertia term (Kh) was adjusted until the modeled 
response curve obtained the correct shape. Adjustment was then made to the lagging 
time constant and gain until the modeled motor power obtained a steady state value 
corresponding to that observed during the response testing. Final adjustment was made 
to the time constant to provide the response shown in Figure 21. 

The final value for the modeled motor power is near but not equal to that observed. 
Some of this variation results from the linear estimation of the mechanical power loss. 
The ideal power loss due to rotational damping is given as [10]: 

Pi^oss = Bat (11) 
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Utilizing per unit values and considering small speed fluctuations near a value of 1, the 
power loss was assumed to behave in a linear fashion with the constant B approximated 
as the average of the two power losses observed during the up-power and down-power 
response testing. 


Open-Loop Response to Increased Electric Load 



Figure 21: Uncontrolled Response to Increased Loading 


Utilizing values identified from the up-power response, the model shown in Figure 19 
was utilized to describe the speed response to a decrease in electric loading. The one 
notable difference between these two events is the sign of the prime mover gain (Kpm). 
For a sudden drop in loading, the increasing motor speed results in a greater counter 
voltage and subsequent drop in armature current. The prime mover gain is negative in 
this case to account for the decrease in power that occurs. Additionally, some adjustment 
to the magnitude of the prime mover gain was required to obtain the response shown in 
Figure 22. Values utilized in the open-loop response models are provided in Table 5. 
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Open-Loop Response to Decreased Electric Load 



Figure 22: Uncontrolled Response to Decreased Loading 

The open-loop response to a change in electric loading represents the self regulating 
aspect of the DC motor as discussed in Section 2.3 where the changing speed affects 
counter voltage and subsequently armature current. Adjustment of terminal voltage is 
required in order to maintain motor speed under varying loading conditions. From 
Equations 5 and 6, the power of the DC motor may be written as: 


M 


\ K 


1 

to 

o 

_1 

D ^ Term ® 

_ M 


— or 

_ R m 


( 12 ) 


The second term of Equation 12 is accounted for in the open-loop response model to a 
change in loading shown in Figure 19. The first term is included to model the 
uncontrolled response to a change in terminal voltage as shown in Figure 23. 
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Vterm Gain Vterm (volts) 


Figure 23: Open-Loop Response Model to a Change in Terminal Voltage 


The affect that a change in terminal voltage has on motor power is modeled directly as 
the product of terminal voltage, constant (Ktv), and speed. Units were introduced in this 
area of the model to assist with proper identification of the closed-loop controller gains 
developed in Section 4.3. The initial estimate for the terminal voltage gain (Ktv) was 
taken directly from Equation 12 as the motor constant divided by armature resistance. 


The power loss gain (Km) was calculated for this model as that loss actually observed at 
the beginning of the test run. Remaining parameters were assigned preliminary values 
based on those used to model the up-power transient since both of these response tests 
have the same initial electrical loading. Adjustment was then made to the prime mover 
gain (Kpm) to maintain the shape of the response curve while the terminal voltage gain 
(Ktv) was manipulated to obtain the correct change in speed. The resulting response 
curve is shown in Figure 24. 
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Open-Loop Response to Increased Terminal Voltage 



Figure 24: Uncontrolled Response to Change in Terminal Voltage 


Table 5 shows parameter values utilized in the three open-loop response models. The 
power loss gain for the step change in terminal voltage was significantly different from 
that observed for step changes in power. This may be a result of variation in the 
armature's interaction with the field given the generator output voltage was initially 118 
volts line to neutral vice the 120 volts established at the beginning of the change in power 
experiments. It is expected that the power loss term will remain near constant with speed 
and voltage control in place. 
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Parameter 

Increased Electric 
Load 

Decreased Electric 
Load 

Increased Terminal 
Voltage 

Kh 

3.3 

3.3 

3.3 

Km 

0.1062 

0.1062 

0.0499 

Kpm 

0.0872 

-0.0805 

-0.0741 

Tpm 

0.144 

0.144 

0.144 

Ktv 

- 

- 

0.262 


Table 5: Open-Loop Model Parameters 


The power loss term (Km) is not actually needed for the model to describe speed 
response; it was included to provide insight into the behavior of modeled motor power to 
assist with the empirical modeling. In lieu of estimating power loss and observing motor 
power, the model could have been constructed utilizing an initial mechanical input equal 
to the initial electrical power. Additionally, as mentioned in Section 3.3, rotational 
friction and windage are often ignored in modeling since they have little influence on the 
performance of the control system. As such, variation in power loss is not considered a 
hindrance to establishing speed and voltage control. 


As shown in Table 5, the inertia term (Kh) and prime mover time constant (Tpm) were 
found to be consistent throughout the three models. The much larger value for Kh, when 
compared to Gh used in the GTG model, indicates that either the machine maintains very 
little energy in an inertial capacitive form as discussed in Section 3.3 or some other 
system dynamic is included in this parameter. Given the large amount of energy lost to 
rotate the generator while maintaining voltage, some of this discrepancy may attest to the 
general imbalance of the tabletop generator. Variation may also be attributed to the 
inherent differences between the DC motor and gas turbine prime movers. For example, 
removal of power from the operational DC motor results in the inertia of the generator 
acting to rotate the motor and in turn causes the motor to behave as a generator with an 
associated counter torque. 
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The variation in the prime mover gain (Kpm) required to fit all three response curves 
presents some concern for the modeling and subsequent implementation of closed-loop 
control. This is addressed in Section 4.3 by evaluating how deviation in this term affects 
the response of the closed-loop system. 

4.3 Closed-Loop Model of Tabletop Generator 

Incorporating what was learned from the open-loop response modeling, closed-loop 
feedback control of the tabletop generator was included and adjustment made such that 
the closed-loop response approximates that of the GTG model developed in Section 4.1. 
The model, shown in Figure 25, follows directly from the open-loop model shown in 
Figure 23 with the addition of a speed error signal input to a PI controller that adjusts 
terminal voltage to drive the speed error to zero. 

Units were included in the modeled controller to assist with exporting gains to the 
software controller. As developed, the product of the frequency error and controller gain 
results in a fractional change in terminal voltage. This fractional change, multiplied by 
the voltage range of the power supply, identifies the desired change in terminal voltage. 

The speed derivative is evaluated in the determination of the sign of the prime mover 
gain (Kpm). As discussed in Section 4.2, this parameter is positive when speed is 
decreasing and negative with increasing speed. This introduces a discontinuity such that 
the default Ordinary Differential Equation (ODE) solver would not converge for a wide 
range of controller gains. In lieu of Simulink's default variable step Domand-Prince 
(ODE-45) solver, the fixed step Runge-Kutta (ODE-4) solver was used with a step size of 
0.001 and a run time of 5 seconds during adjustment of controller gains. 

Utilizing parameters indentified in Table 5 for an increase in electric loading from 0.045 
to 0.135 per unit and the terminal voltage gain (Ktv), controller gains were adjusted until 
the response resembled that of the GTG model as shown in Figure 26. After reaching the 
desired response, the terminal voltage and controller gains were adjusted correspondingly 
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to provide a fractional change in terminal voltage matching that observed during similar 
response testing. Values utilized in this model are shown in Table 6. 



Fractional 
Change in Vterm 


Figure 25: Tabletop Closed-Loop Response Model 
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Parameter 

Value 

Kp 

0.165 

Ki 

0.215 

Kh 

3.3 

Km 

0.1062 

Kpm 

0.0872 

Tpm 

0.144 

Ktv 

0.36 

Cl 

274.2 

C2 

150 

C3 

0.0002 


Table 6: Tabletop Closed-Loop Parameters 


Response to Electric Load Step Change (0.045 to 0.135) 



Time (sec) 


Figure 26: Tabletop Model Closed-Loop Response 
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Given the uncertainty in the estimation of the power loss and the variation in the prime 
mover gain seen in open-loop modeling, the closed-loop model was revised as shown in 
Figure 27. With estimation of power loss removed from the model, the initial steady 
state condition was established by including an initial mechanical power such that the 
power mismatch is zero at time equal to zero. Using the same 0.045 to 0.135 step 
change, the prime mover gain was adjusted upward and downward 15% to observe how 
sensitive the response is to deviations in this parameter. Results of this analysis are 
shown in Figure 28. 



volts 


Gain 


Figure 27: Tabletop Closed-Loop Model without Estimation of Power Loss 
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Comparison of the response curve for the non-deviated prime mover gain (1.0 Kpm) 
shown in Figure 28 with that shown in Figure 26 provides validation that rotational 
energy losses (windage, friction, etc.) need not be included in the closed-loop model 
since they have little affect on the controlled response. The slight variation in modeled 
response, given a 15% deviation in the prime mover gain, attests to the stability and 
robustness of the closed-loop feedback control. 


Tabletop Model Sensitivity to Prime Mover Gain 



Time (sec) 


Figure 28: Tabletop Model Sensitivity to Prime Mover Gain 


The tabletop model evaluation to this point has been focused on the response to a specific 
change in electrical loading of 0.045 to 0.135 per unit. This was set as the design point to 
correspond to the readily available 225 to 675 watt step loading of the tabletop generator 
rated at 5 KW. Using the Simulink model with no estimation of power loss, the response 
to a step change in electric loading from 0 to 0.225 is shown in Figure 29. This step in 
power represents the uppermost power excursion for the tabletop generator which is 
limited in its useable power output to a little more than 1125 watts. 
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Response to Electric Load Step Change (0 to 0.225) 



Time (sec) 


Figure 29: Closed-Loop Response to a Larger Step Input 

As shown, utilizing controller gains identified for the step change in power from 0.045 to 
0.135 per unit does not provide the same level of agreement when a larger step change of 
0.225 per unit is applied to the gas turbine and tabletop models. Adjustment of the 
tabletop controller's proportional gain from 0.165 to 0.19 improves the correlation 
between peak values. 
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Chapter 5 


Hardware Model 

A physical description of the hardware model is provided in Chapter 5 along with 
response testing results. Appendix D provides a listing of parts, including specifications 
and purchasing information. Recommendations to improve the response of the tabletop 
generator are made based on analysis of actual speed sensor noise. 

5.1 Physical Description 

The hardware model of a shipboard generator as constructed is pictured in Figure 30. 

The 3-phase 5 KW generator is side driven by a timing belt with a 32 tooth pulley 
mounted on the generator shaft and a 22 tooth pulley on the motor. Two DC motors, in 
series electrically and mechanically coupled, provide a rated capacity of 4 HP to serve as 
the prime mover. Not seen in the picture is a small cooling fan that provides forced 
circulation to the motors and a Hall Effect speed sensor rigidly mounted adjacent to a 14 
tooth sprocket set on the furthermost motor shaft. Two XHR 1000 Watt Series 
programmable DC power supplies provide a rated 14 amps and 150 volts to the DC 
motors. The additional power supply pictured provides manual control of current to the 
generator's field windings in order to maintain an output voltage of 120 V line to neutral 
(208 V across phases). A programmable power supply for voltage control was not in use 
at the time of this photo. 
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Figure 30: Tabletop Generator as Constructed 


The analog / digital USB interface board used to pass measured and controlling signals to 
and from the software environment is shown in Figure 31. Additional circuit boards were 
constructed as necessary to modify speed and voltage measurements prior to the USB 
interface. Labview was utilized to implement speed and voltage control. 



Figure 31: USB Interface Board 
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Implementation of the control scheme in the software environment was conducted by 
Vanessa Esch while working on an Undergraduate Advanced Project. Signal 
conditioning of speed and voltage measurements was developed by Jacob Osterberg 
while working through the Undergraduate Research Opportunity Program. Their work is 
shown in its entirety in the Appendices. Though some information is duplicated and 
conflicts may exist, no attempt was made to modify or summarize their reports. This 
preserves the insight of three observers involved on one project. 

The loading scheme utilized for response testing is shown in Figure 32. A bank of 9 light 
fixtures allows for some number of variations in electric loading. Though the resistance 
of a light bulb changes to some degree with temperature, the purely resistive loading was 
preferred for its ease and reasonable representation of step changes in power. The 
fixtures are wired in a Y configuration with each light seeing a voltage of 120 V line to 
neutral. Each column of fixtures is supplied power from one phase. 

The first row of lights is operated via the three rightmost switches. Operation of these 
switches, as near to simultaneous as possible to minimize unbalanced loading, allows for 
establishment of an initial loading condition if desired. The second and third rows are 
provided power via the three relays pictured in the left corner. Operation of the leftmost 
switch energizes the relays to provide instantaneous loading in order to model a step 
change in power. 
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Figure 32: Loading Scheme for Response Testing 


General characteristics of the tabletop generator are summarized in Table 7. The 
inclusion of a saturation block in the closed-loop modeling of Section 4.3 was considered 
but deemed impractical. The saturation represents the limitation of the power supply and 
was initially thought to be a simple limit placed on the terminal voltage in the model. As 
observed, a sudden increase in terminal voltage, as required to maintain speed for a large 
up-power transient, drives the power supply to a current limited condition for some 
period of time. This occurs while increasing power from 225 watts to some value near 
750 watts. As such, there is some limit on up-power transients that may be conducted 
while maintaining good speed response behavior. 
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Parallel Power Supply Rating: 

2100 watts (14A, 150V) 

Power Supply Usable: 

2025 watts (current limited) 

Motor Speed: 

2618 rpm 

Motor Constant: 

0.432 volt*sec (motors coupled) 

Motor Resistance: 

1.87 ohms (motors coupled) 

Power Lost to Heat: 

365 watts 

Motor Power: 

1660 watts 

Generator Speed: 

1800 rpm 

Power Loss: 

~ 525 watts 

Generator Output: 

-1135 watts (tested up to 1125 watts) 


Table 7: General Characteristics of Tabletop Generator 


The power loss of the generator is significant. About half of this value is attributed to 
windage and rotational friction. The remaining loss is observed upon energizing the field 
windings from an external source to provide required output voltage. Utilization of the 
generator's installed voltage regulator vice an external power source exacerbates the 
power loss. 

5.2 Response Testing 

Response testing proved to be very time consuming and did not yield desired results. 
Electromagnetic interference introduced significant fluctuations in the speed sensor 
signal. This occurred after the generator was electrically loaded for a short period of 
time. Relocation and shielding of the speed sensor wiring alleviated this problem. The 
best response of the tabletop generator is not similar to the GTG model as shown in 
Figure 33. 
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Tabletop Generator Response 



Time (sec) 

Figure 33: Response of Tabletop Generator 

The proportional and integral gains in the software controller were set to 0.0012. Values 
higher than this resulted in significant oscillation. As shown previously in Table 6, 
modeled proportional and integral gains to match the GTG response are 0.165 and 0.215. 
The proportional and integral gains of the tabletop model were set at 0.012 to match the 
experimental response as shown in Figure 33. Given the similar response obtained with 
gains 10 times greater than those actually used, it is possible that a factor of 10 error 
exists between the Simulink model and the Labview controller. 

Speed signal noise is a likely contributor in the significant oscillation that occurs with 
increased controller gains. The filtered speed signal is compared to raw data in Figure 34 
and Figure 35. The effect of the signal's 8-bit resolution is shown in Figure 36. 
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Filtered Speed Signal 



Time (sec) 


Figure 34: Filtered Speed Signal 


Filtered Speed Signal 



Figure 35: Magnified View of Filtered Speed Signal 
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8-bit Affect on Speed Signal 
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Figure 36: 8-bit Resolution of Speed Signal 

The 8-bit resolution of the speed signal results in a non-uniform oscillation of the speed 
error where one value tends to be favored. As shown above, the signal is artificially low. 
The amplitude of the filtered speed noise is slightly greater than the 8-bit signal but 
shows more centered values. The 8-bit signal noise was introduced to the tabletop model 
as a repeating sequence similar to that shown in Figure 36 where one value is favored 
over another. The modeled response is shown in Figure 37 with the required controller 
gains to match the GTG model. 



The oscillations predicted by the model occur at a rate faster than the speed sensor data 
rate. Currently, 14 data points are averaged to represent each motor revolution. This 
results in approximately 43 signals per second or one signal every 0.023 seconds. The 
steep oscillations predicted by the model suggest a significant change in speed may occur 
over this duration. The modeled affect of increasing controller gains is shown in Figure 
38 and Figure 39. 
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Model with 8-bit Resolution Signal Noise 



Time (sec) 


Figure 37: Tabletop Model with 8-bit Signal Noise 


Increasing Tabletop Model Gains with 8-bit Signal Noise 



Figure 38: Increasing Controller Gains 




























































Increasing Tabletop Model Gains with 8-bit Signal Noise 



Figure 39: Increasing Controller Gains, Magnified 


As gains are increased, the modeled response to speed signal noise shows steep 
oscillations that may not be seen by the controller until a significant change in speed has 
occurred. As such, speed sensor noise must be reduced while increasing speed data rate. 
The 14 tooth sprocket currently installed is not designed specifically to interface with the 
Hall Effect speed sensor. Construction of a 'speed sensor sprocket' may reduce noise and 
increase data rate. 

Increased resolution of the speed signal may also improve response. However, it is 
possible that the 8-bit resolution may provide acceptable performance if signal noise is 
adequately reduced. Inclusion of gain scheduling is another option for consideration. 
Sensing a change in electrical loading, larger controller gains would be utilized during the 
transient while smaller gains are used during steady state operation. 
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Chapter 6 


Conclusions and Future Work 

The numerical model developed to describe the response of the Allison 501-K34 gas 
turbine generator provides the key characteristics of the actual response. It accurately 
describes the peak values reached during a transient and the time required to reach steady 
state. It does not precisely model the transition from the peak value to steady state. 

The closed-loop feedback control modeling of the tabletop generator strongly suggests 
that the tabletop generator driven by a DC motor may be controlled in a fashion to 
resemble the speed response of a gas turbine generator. The best agreement occurs when 
the controller gains are adjusted for a specific change in loading. Deviations from this 
step change in electric power results in a modeled response of the tabletop generator that 
reaches a peak value slightly different than that predicted by the gas turbine model. 

The response of the tabletop generator as realized is not similar to the GTG model. 
Controller gains utilized are significantly lower than those suggested by the modeling to 
match the GTG. As gains are increased, the modeled response to speed signal noise 
shows steep oscillations that may not be seen by the controller until a significant change 
in speed has occurred. 

The model gains required to match the actual response of the tabletop generator are off 
by one order of magnitude. A factor of 10 error may exist between the model and 
software controller. 

Speed sensor noise must be reduced while increasing speed data rate. The 14 tooth 
sprocket currently installed is not designed specifically to interface with the Hall Effect 
speed sensor. Construction of a 'speed sensor sprocket' may reduce noise and increase 
data rate. 
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Increased resolution of the speed signal may also improve response. However, it is 
possible that the 8-bit resolution may provide acceptable performance if signal noise is 
adequately reduced. 

Gain scheduling is another option for consideration. Sensing a change in electrical 
loading, larger controller gains would be utilized during the transient while smaller gains 
are used during steady state operation. 

Future work for this project includes that necessary to model a shipboard electrical 
distribution system. This includes developing a control system that provides accurate 
load sharing amongst two or more isochronous generators and establishing a 
methodology to parallel these generators in phase. 

Additional work may also include the incorporation of an angled Howse Gearbox that has 
a well suited gear ratio for this application. This may reduce some power lost to 
rotational friction associated with the side loading of the generator shaft by the timing 
belt drive. 
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Appendix A: Nomenclature, Data, and Calculations 

1) Nomenclature as used in this thesis is listed below. 


Cl 

conversion constant, per unit speed to rad/sec 

C2 

conversion constant, fractional change in Vrerm to volts 

C3 

conversion constant, watts to per unit power 

Er 

error 

fmax 

maximum fuel flow 

fmin 

minimum fuel flow 

Gh 

inertia term, GTG model 

H 

inertia constant 

Im 

motor current 

K 

motor constant 

Kc 

controller gain 

Kdg 

governor droop 

Kgl 

governor gain 

Kg2 

governor gain 

Kg3 

governor gain 

Kh 

inertia term, tabletop model 

Ki 

integral control gain 

Km 

mechanical power loss gain 

Kp 

proportional control gain 

Kpm 

prime mover gain 

Ktv 

terminal voltage gain 

OP 

controller output 

Pm 

motor power 

Rm 

motor resistance 

Td 

derivative time constant 

Tfl 

fuel valve time constant 

Tgl 

governor time constant 

Tg2 

governor time constant 

Tg3 

governor time constant 

Tg4 

governor time constant 

Ti 

integral time 

Tm 

motor torque 

Tpm 

prime mover time constant 

Ttl 

turbine leading time constant 

Tt2 

turbine lagging time constant 

V c 

counter voltage 

Vxerm 

terminal voltage 

CO 

rotational speed (rad/sec) 
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2) Data used in the determination of motor constants. Counter voltage was plotted as a 
function of speed (rad/sec) where the slope of the resulting line equals the motor constant 
K (volt*sec) as discussed in Section 2.3. 


Voltage (V) 

Current (A) 

Motor 2 Driving Motor 1 

Speed (rpm) Counter Voltage (V) 

Speed (rad/sec) 

10 

1.14 

382 

8.6 

40.0 

20 

1.21 

824 

18.6 

86.3 

30 

1.29 

1265 

28.5 

132.5 

40 

1.34 

1714 

38.7 

179.5 

50 

1.38 

2152 

48.6 

225.4 

60 

1.42 

2599 

58.7 

272.2 

70 

1.48 

3044 

68.7 

318.8 

80 

1.53 

3489 

78.6 

365.4 

90 

1.54 

3932 

88.5 

411.8 

100 

1.55 

4380 

98.4 

458.7 

110 

1.59 

4822 

108.2 

505.0 


Voltage (V) 

Current (A) 

Motor 1 Driving Motor 2 

Speed (rpm) Counter Voltage (V) 

Speed (rad/sec) 

10 

1.22 

371 

8.6 

38.9 

20 

1.34 

814 

18.6 

85.2 

30 

1.42 

1250 

28.5 

130.9 

40 

1.48 

1691 

38.6 

177.1 

50 

1.53 

2137 

48.8 

223.8 

60 

1.58 

2579 

58.9 

270.1 

70 

1.64 

3026 

68.9 

316.9 

80 

1.70 

3467 

79.0 

363.1 

90 

1.70 

3919 

89.2 

410.4 

100 

1.74 

4366 

99.3 

457.2 

110 

1.69 

4817 

109.5 

504.4 


Resulting Motor Constants: Ki = 0.215 (volt*sec) 

K 2 = 0.217 (volt*sec) 
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3) Terminal voltage and speed data collected from both motors operating in series with 
no mechanical loading. Note that a 200 volt power supply was initially used for some 
testing prior to utilizing the 150 volt programmable power supplies. 


Motor 1 and 2 in Series, No Load 
Voltage (V) Current (A) Speed (rpm) Speed (rad/sec) 


20 

0.57 

418 

43.8 

40 

0.59 

862 

90.3 

60 

0.62 

1309 

137.1 

80 

0.65 

1750 

183.3 

100 

0.67 

2200 

230.4 

120 

0.70 

2650 

277.5 

140 

0.72 

3096 

324.2 

160 

0.73 

3546 

371.3 

180 

0.74 

3995 

418.4 

190 

0.77 

4226 

442.5 


4) Motor Resistance 



Motor 2 Driving Motor 1 



Voltage to „ , , . 

Drive Motor (V) Speed < rpm > 

Driven Motor 

Driven Motor 


Counter Voltage (V) 

Current (A) 

Driven Motor Open 

5 166 

3.79 

- 

Driven Motor Shorted 

8.1 166 

- 

3.96 



Motor 1 Driving Motor 2 



Voltage to 

Speed (rpm) 

Driven Motor 

Driven Motor 


Drive Motor (V) 

Counter Voltage (V) 

Current (A) 

Driven Motor Open 

5 

166 

3.83 

- 

Driven Motor Shorted 

8.6 

166 

- 

4.2 


3.79 
•Mi - 3 96 

R M1 = 0-957 

3.83 
« '* 4.2 

R M2 = °-912 
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5) Torque-speed curve for both motors operating in series electrically while 
mechanically joined via a rigid coupling. 



Speed (rpm) 

Vterm (volts): 

Speed (rad/sec) 

100 

Torque (N*m) 

125 

Torque (N*m) 

150 

Torque (N*m) 

K (V*sec) 

0 

0.0 

23.1 

28.9 

34.7 

0.432 

400 

41.9 

18.9 

24.7 

30.5 


800 

83.8 

14.7 

20.5 

26.3 

Rm (ohms) 

1200 

125.7 

10.6 

16.3 

22.1 

1.87 

1600 

167.6 

6.4 

12.2 

17.9 


2000 

209.4 

2.2 

8.0 

13.8 


2400 

251.3 

-2.0 

3.8 

9.6 


2800 

293.2 


-0.4 

5.4 


3200 

335.1 



1.2 


3600 

377.0 



-3.0 


Sample Calculation: 


K := 0.432 (volfsec) 
R M = 1 .87 (ohms) 
v Term ;= 100 ( volts ) 


o := 


f 0 

41.9 

83.8 

125.7 

167.6 

209.4 

Ul-3y 


(rad sec) 



Term 


|V N 

R M. 


CO 


(N*m) 
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6) Best Suited Motor Speed and Gear Ratio 


V Term := 150 (volts) K := 0 432 ( volt * sec ) 
I M = 14(amps) R M := 1.87 (ohms) 


('60') 

Speed := co — Speed = 2737 (rpm) 

\2-V 


co := 


'’Term “ ! M R M 


K 


co = 286.6 (rad sec) 


2737 

Gear Ratio := - Gear Ratio =1.52 

1800 


7) Analysis of Available Gear Ratios 
Power Supply Voltage Limited 



' 1.6 " 




^301.6 > 

Gear_Ratio := 

1.5454 

co := Gear_Ratio 1800 


co = 

291.3 


, 1-52 , 


^ 60 J 


^286.5 , 


(rad sec) 


v Term :=150 ( volts ) I\f = 


M 


v Term “ K “ 
R M 


( 10.5' 

I M = 12-9 
l H . 


(amps) 


T 1581' 

PwrSup := I^|-V>p erm PwrSup = 1938 

Ul04; 


(watts) 


Power Supply Current Limited 



r 1.52 ' 




( 286.5'' 

Gear_Ratio := 

1.47 

co := Gear_Ratio 1800 


co = 

277.1 


^ 1.4545 , 


s 60 J 


G>74.2, 


(rad sec) 


I]VI ;= I 4 (^P 5 ) 


v Term = ! M R M + Kco 


( 150 ' 
VTerm = 145-9 
U44.6. 


(volts) 


Pwr_Sup := I M V Term 


f 2099^ 

Pwr_Sup = 2042 

12025 , 


(watts) 
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8) Determination of the Inertia Constant (H) and the inverted inertia term (Gh) utilized in 
the Simulink model of the Allison 501-K34 GTG. 

The inertia constant is defined as the Kinetic Energy of the rotating system divided by the 
Volt-Amp rating of the generator. It represents the amount of time that the generator may 
supply power at its full rating utilizing the capacitive nature of the inertia with no 
additional mechanical input. 

Note that the shipboard generator is a 4 pole machine that operates at 1800 rpm and it is 
not known if the inertia of the generator includes the associated gearbox. 


I GT := 61 (lb*fE) 1 ^ := 3550 (lb*ft 2 ) I sys := (l GT + 1 ^) (0.04214) 


I svs = 152.2 (kg*m 2 ) o := 1800- 


' 2-0 


V 60 y 


co = 188.5 (rad sec) 


1 2 6 
KE S y S := - Tgyg Co VA_Rating := 3 10 (volt amps) 


KE S y S 

VA_Rating 


H = 0.9 (sec) 



Gh = 0.555 (sec 4 ) 


9) Tabletop Open-Loop Response Testing 

Note that at the time of this testing, the prime mover gear ratio was less than optimum 
with a 32 tooth pulley on the generator and a 20 tooth pulley on the motor resulting in a 
gear ratio of 1.6 and a motor speed of 301.6 rad/sec to provide 60 Hz. As discussed at the 
end of Section 2.3, a better choice of a 22 tooth pulley was later installed on the motor 
resulting in the finalized motor speed of 274.2 rad/sec utilized elsewhere. Determination 
of initial motor power and power loss follows. 

Supplying 225 watts at 60 hz and 12Q volts (line to neutral) 

K := 0.432 (volt*sec) I M := 5.85 (amps) co := 1.0056 (301.6) (rad/sec) 


P\f := K lyj co Pyj = 766 (watts) ^loss := P M — 225 Pj oss = 541 (watts) 

P M 

P M := 5000 p M = 0-1533 (per unit) 
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Supplying 675 watts at 60 hz and 120 volts (line to neutral) 


K := 0.432 (volt*sec) I M := 9.21 (amps) a := 0.9991 (301.6) (rad sec) 

P M ;=KI M® P M = 1199 (watts) p loss := P M“ 675 P loss = 524 ( watts ) 

P M 

P M : =™ p M = 0.2398 (per unit) 


Supplying 225 watts at 58.2 hz and 118 volts (line to neutral) 

K := 0.432 (volt*sec) I M := 3.74 (amps) a := 0.9738 (301.6) (radsec) 

p M := K co P M = 475 (watts) P loss := P M - 225 P loss = 250 (watts) 

P M 

P M := 5000 P M = 0 0949 (per UfUt) 

The raw speed data required filtering to provide a reasonable response curve. Averaging 
the forward 10 data points in Excel resulted in the following improvement. 
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Sample data from open-loop response testing is provided below. The actual measured 
parameter is the frequency of the square wave generated as a sprocket mounted on the 
motor shaft interacts with the speed sensor fixed in close proximity. The inverse of the 
measured frequency is the period or pulse duration. Summation of the pulse durations 
provide the indicated time. With 14 teeth on the sprocket, the following converts sensor 
frequency to motor speed: 



Speed 

= Freq ■ 

( 2^ 

V 14 ; 


Speed Sensor Frequency 

Pulse Time 

Time 

Motor Speed (rad/sec) 

Filter (fwd avg 10) 

686.44 

0.00146 

0.5521 

308.0718 

303.93 

669.16 

0.00149 

0.5536 

300.3205 

303.26 

671.32 

0.00149 

0.5551 

301.2882 

303.41 

679.53 

0.00147 

0.5566 

304.9735 

303.36 

676.04 

0.00148 

0.5581 

303.4065 

303.25 

680.64 

0.00147 

0.5595 

305.4717 

303.46 

683.06 

0.00146 

0.5610 

306.5567 

303.25 

676.04 

0.00148 

0.5625 

303.4065 

302.58 

671.50 

0.00149 

0.5640 

301.3692 

303.10 

678.24 

0.00147 

0.5655 

304.3943 

303.10 

671.50 

0.00149 

0.5669 

301.3692 

303.21 

672.59 

0.00149 

0.5684 

301.8556 

303.62 

670.24 

0.00149 

0.5699 

300.8036 

303.46 

677.14 

0.00148 

0.5714 

303.8996 

303.72 

680.64 

0.00147 

0.5729 

305.4717 

303.45 

676.04 

0.00148 

0.5743 

303.4065 

302.64 

668.09 

0.00150 

0.5758 

299.8390 

302.38 

687.57 

0.00145 

0.5773 

308.5801 

302.23 

671.50 

0.00149 

0.5788 

301.3692 

301.20 

680.64 

0.00147 

0.5803 

305.4717 

301.19 

680.64 

0.00147 

0.5817 

305.4717 

300.48 

669.16 

0.00149 

0.5832 

300.3205 

299.66 

675.86 

0.00148 

0.5847 

303.3245 

299.71 

671.32 

0.00149 

0.5862 

301.2882 

299.61 

662.43 

0.00151 

0.5877 

297.2966 

299.61 

670.24 

0.00149 

0.5892 

300.8036 

299.86 

664.72 

0.00150 

0.5907 

298.3242 

299.51 

664.72 

0.00150 

0.5922 

298.3242 

299.96 

671.32 

0.00149 

0.5937 

301.2882 

299.86 

664.72 

0.00150 

0.5952 

298.3242 

299.46 

662.43 

0.00151 

0.5967 

297.2966 

299.06 

670.24 

0.00149 

0.5982 

300.8036 

299.41 

673.67 

0.00148 

0.5997 

302.3437 

298.97 

671.32 

0.00149 

0.6012 

301.2882 

298.77 


Sample Data from Increased Electric Loading 
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Data collected prior to and after the desired response was removed after close 
examination of the areas of interest. The following shows how the data was shifted from 
when the data was collected (time stamp) to the actual response time. Also shown is the 
determination of initial speed which comes from the averaging of data points prior to the 
response time of zero. 




Fwd Avg 10 



Time Stamp 

Response Time (sec) 

Speed (rad/sec) 

Speed (per unit) 



0.5714 

0.0000 

303.72 

1.0070 

Initial Speed 

0.5729 

0.0015 

303.45 

1.0061 

303.28 

1.0056 

0.5743 

0.0029 

302.64 

1.0034 

303.05 

1.0048 

0.5758 

0.0044 

302.38 

1.0026 

304.30 

1.0089 

0.5773 

0.0059 

302.23 

1.0021 

303.63 

1.0067 

0.5788 

0.0074 

301.20 

0.9987 

303.78 

1.0072 

0.5803 

0.0089 

301.19 

0.9986 

304.23 

1.0087 

0.5817 

0.0103 

300.48 

0.9963 

303.71 

1.0070 

0.5832 

0.0118 

299.66 

0.9936 

303.49 

1.0063 

0.5847 

0.0133 

299.71 

0.9937 

302.82 

1.0041 

0.5862 

0.0148 

299.61 

0.9934 

302.44 

1.0028 

0.5877 

0.0163 

299.61 

0.9934 

302.74 

1.0038 

0.5892 

0.0178 

299.86 

0.9942 

303.03 

1.0047 

0.5907 

0.0193 

299.51 

0.9931 

302.38 

1.0026 

0.5922 

0.0208 

299.96 

0.9946 


1.0056 

0.5937 

0.0223 

299.86 

0.9942 



0.5952 

0.0238 

299.46 

0.9929 



0.5967 

0.0253 

299.06 

0.9916 



0.5982 

0.0268 

299.41 

0.9928 



0.5997 

0.0283 

298.97 

0.9913 



0.6012 

0.0298 

298.77 

0.9906 



0.6027 

0.0313 

298.22 

0.9888 



0.6042 

0.0328 

297.77 

0.9873 



0.6057 

0.0343 

298.37 

0.9893 




Sample of Response Data Utilized in Open-Loop Modeling 
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10) Employed Speed Signal Filtering 

The speed signal as utilized employed a filter that averaged every 14 data points. This 
corresponds with one speed signal for each motor rotation. The Matlab code used to 
create Figures 34, 35, and 36 is shown below. 


UF = dlmread('Unfiltered_SS.txt'); %UnFiltered data at Steady-State 
%column 1 is time, column 2 is per unit speed 

for i=l:(length(UF)/14) 

Avgl4(i,l) = UF(1+14*(i-1),1); %time at every 14th data point 
for j=l:14 

speed(j) = UF(j+14*(i-1),2); 

Avgl4(i,2) = mean(speed); %avg every 14 speed values 

end 

end 

%discrete values for speed 
%range is 50-70hz, 8-bit 

Res = (20/256)/60; %resolution of per unit speed 
for i=l:(length(UF)/14) 

Dis(i,l) = Avgl4(i,l); %time 

Dis(i,2) = round(Avgl4(i, 2)/Res)*Res; %discrete speed value 

end 

hold on 

plot(UF(:,1),UF(:,2)) 

plot(Avgl4(:,1),Avgl4(:,2), 'k' ) 

%plot(Dis(:,1),Dis(:,2)) 
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Appendix B: Implementation of the Control System, by Vanessa Esch 


Control System for a DC Motor-Generator System 
for operation as a Shipboard Gas Turbine Generator Model 


Vanessa Esch 
6.UAP 
Spring 2009 
at 

Laboratory of Electromagnetics and Electronics 
Massachusetts Institute of Technology 


Abstract 

The operation of a U.S. Navy destroyer power system is modeled in this project. The gas 
turbine generator system has been modeled using a DC motor and generator. The goal of the 
project is to prove that this model can be extended to act as a scaled GTG system on a TableTop 
platform to test load faults of the shipboard system. The model was completed first in Simulink 
and then a physical system to control the model was implemented. The work of the author of 
this paper was on implementing the physical control system. Although currently the system has 
power limitations, it can successfully be controlled, and with higher power supplies, can mimic 
the response of the Navy GTG system. 
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1. Introduction 

The purpose of this project is to develop a physical system that operates under reduced 
scaled characteristics of a U.S. guided missile destroyer (DDG) gas turbine generator set (GTG). 
The scaled system that has been developed uses a DC motor and generator instead of the gas 
turbine. The design goal is to use modeling and simulation to find the appropriate control 
constants for the system such that it behaves as a reduced scale DDG plant would. The control 
system is then implemented as a software-based feedback system for operation of the DC motor- 
generator set. Once this scaled system is fully proven and tested via implementation, the next 
step for the project is to do load and fault analysis in a manner similar to what is demanded of the 
system a DDG. The modeling and simulation aspect of the project is the work of a U.S. Naval 
Engineering Duty Officer thesis project. This undergraduate advanced project is to design and 
implement the actual control loop for the system in software. 

This paper presents the technical background and work done toward these goals and the 
particular contribution of the author. Section two provides a technical overview of the GTG 
system on a DDG-51 class ship including its method of feedback control and the operating 
requirements of the shipboard system. Section three provides a similar description of the scaled 
system developed for testing in the project. This overview description will be followed by an in- 
depth analysis of the design decisions made for the system in section four. The section will 
focus on the control system, including the algorithm, physical system, and user interface, as that 
was the primary work of the author of this project. Section five is a complete explanation of how 
to use the system and user interface in the hope that it can serve as the foundation for future 
testing and use of the scaled system. Section six will conclude the paper with a reflection on the 
project and suggestions for future development of the system. 
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2. Navy System 

2.1 Destroyer Physical System 

Power generation of a DDG-51 class ship is provided by three 2500kW gas turbine 
generator (GTG) sets. The gas turbines are LM2500 engines made by General Electric. Each 
GTG set produces three-phase 450 VAC 60Hz for shipboard operations. 1 The Allison 501-K34 
Gas Turbine Engine is the prime mover on the DDG-51 class ship; it is a single shaft, axial flow 
gas turbine. For 60Hz output from the generator, the LM2500 operates at 14,340 RPM. 
Reduction gear translates the engine output to 1800RPM, which correlates with 60Hz output. 

2.2 Control System 

On the highest level, the control system implemented on DDG GTG systems consists of 
two elements, a governor and a voltage regulator. These two components are very commonly 
used in control systems for prime mover and generator sets. 11 The governor is the element that 
regulates the output frequency of the system and the voltage regulator controls the amplitude of 
the voltage produced by the generator. A basic understanding of the design of these elements on 
the naval system directed the design of the system in the scaled model. The description of these 
systems is from training documents provided by the Surface Warfare Officers School in 
Newport, RI. 1 " 

On a generator system, real power is related to the prime mover’s delivered energy, and 
therefore the speed of the generator. Reactive power is related to terminal voltage. A generator 
fails under too great of a real power load due to speed collapse and also fails under too great a 
reactive load due to voltage collapse. lv The system built here attempts to avoid both types of 
collapse under scaled, rated load application and removal. 


80 



3 

2.2.1 Governor 

A change in the ship’s load is sensed by a load sensor at the output of the generator. 

When a change in the load occurs, the output voltage and current signals will change both in 
amplitude and frequency. v In order to regulate the frequency of the system, the load sensor 
measures the frequency of the generator output voltage and calculates a speed error signal for the 
feedback system by comparing that value to a speed reference value set by the operator, 
nominally 60 Hz. vl 

The speed of the prime mover (gas turbine) is regulated by a feedback loop based on this 
speed error signal. The speed error is sent through both stability and gain amplifiers, “the gain 
amplifier decreases response time (dampens); the stability amplifier prevents hunting.” v " In 
basic control theory, Proportional-Integral-Derivative (PID) control is the common type of 
control system, and among the three parallel branches of a PID control, integral control is the 
primary mechanism for avoiding drooping and derivative is the primary mechanism for 
preventing hunting (ie, oscillate) which occurs which pure proportional control.™ 1 A very 
common method of control for GTG found in literature is PI control (without a derivative 
branch).” 1 For the naval system, it is appears that PI control is the means used for closed loop 
feedback in the speed governor. The system used in the model also implements, optionally, a 
dead band to avoid hunting near the set-point. x 

There are two modes of frequency control on the DDG, isochronous and droop (ISO and 
DROOP). xl In ISO operation, feedback will try to maintain frequency at the reference speed 
while in DROOP mode, feedback will not occur and frequency will drop if the load increases. 

On the ship, DROOP mode is used only for connecting to shore power as shore power is 
considered an infinite bus supply and if load sharing between it and the shipboard generators 
takes place, the generators will be pushed beyond their capacity. 
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The steps required to implement the naval speed control system are more complex than 
the ones implemented in the scaled version presented here because of the different prime movers. 
In order to adjust speed, the gas turbine fuel input rate must adjust and temperature bounds must 
be monitored closely. In the scaled version, all that is necessary to adjust the speed of the DC 
motor is an adjustment of the DC input voltage to the motor. 

2.2.2 Voltage Regulator 

The primary purpose of the voltage regulator is to directly adjust the excitation field 
current applied to the rotor of the generator so that the output voltage of the system has the 
desired amplitude.* 11 Additional components in the voltage regulation system on the DDGs 
include a current limiter, an over-voltage limiter, and a failure detection module. 

The voltage regulator is supposed to maintain the following operating restrictions on the 

system; 

1. Current over 12,900A on any phase of the generator stator will result in a 
reduction in field current such that the system voltage (and therefore current) will 
be reduced. 

2. Voltage over 522 VAC for 1.5 seconds or 590 VAC instantaneously will force the 
field current to drop to the level that correlates to 518 VAC output until an over¬ 
voltage protection switch is re-set. 

2.3 Navy Operating Requirements 

The scaled model developed in this project must operate under the constraints of the GTG 
system actually deployed on the DDGs. The naval specifications have strict voltage and 
frequency requirements, these requirements differ slightly between flights (or versions) of the 
DDG class of ships. Flight I ship, DDG 51-90, have different GTG specifications than Flight IIs, 
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DDG 91 and later. For Flight I ships, the application and removal of rated load may cause no 
greater momentary deviation from system frequency than two percent. However, on the Flight 
I's, this maximum error is two percent for overspeed (application of load) and 2.5 for underspeed 
(removal of load).*" 1 

On both flights and for load removal and application, the system has 1.5 second to return 
to within one percent of the normal operating frequency. This final requirement is interpreted to 
mean that during normal operation, frequency deviation within one percent is acceptable. Also, 
any added or removed load that does not cause the steady-state load to exceed the rated load 
should not trip the engine. In order to meet these specifications, a well calibrated control system 
will be required of the scaled system. 

3. Scaled Test System 

3.1 Physical System Design 

3.1.1 DC Motor 

The system controlled in this project consists of two series permanent magnet DC motors 
that drive a generator. The motors are Leeson 2 HP machines; in series they deliver a total rated 
capacity of 4 HP. Basic Kirchoff s laws establish the characteristic equation for a DC motor, see 
equation [1]. The equivalent circuit model of a DC motor is shown in Figure 1. 

V term ^a "l" ^emf (1) 
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Figure 1: Equivalent circuit of DC motor 



In a permanent magnet motor, the armature current is held constant, therefore there is a 
linear relationship between speed and back electro-motive force (EMF) (2).[Greg's 1] Torque 
also exhibits a linear relationship, albeit with armature current (3). Experimentation in the lab 
determined the motor constant, K, to be .215 Volt-sec for one motor and .217 Volt-sec for the 
other. xiv 



(2) 

T = K*l a 

(3) 


An increase in mechanical load will decrease counter voltage. Since terminal voltage can 
remain constant under a change of load, this change would force the armature current to become 
much larger in order to satisfy the 
Kirchoff s voltage loop. The DC 
motor shows its inherit stability 
in this situation because armature 
current and counter voltage 
would adjust until, at a lower 
speed, power could be matched 
between the motor and load. 


Figure 2: Torque-Speed curve ot scaled system 


Torque-Speed for 2 Motors in Series 



84 


















7 

This would be a new steady state for the motor. A change in terminal voltage would have to take 
place in order to get the speed of the system back to its original level. 

The torque-speed curve is the essential operating curve for these operations. Although 
the two motors in series have a rated speed of 4800 RPM at 240 V, the system limitation is 150 
V for terminal voltage. Therefore, torque-speed curve, measured by Greg Elkins, shows terminal 
voltages of 150V or less. The torque-speed curve is shown in Figure 2. 

The generator is a Ming Dong STC-5 series three-phase AC synchronous generator. xv 
The generator has a rotary field and is three-phase Y-connected with a center ground. The 
generator is a salient pole machine. It has line-to-line and line-to-phase voltages of 400 V and 
230 V, respectively. Its rated output is 6.3 KVA or 5KW with a rated current of 9A and power 
factor of .8 lagging. Rated speed is 1500 r/min. The generator has 4 poles. The prime movers 
for these generators can be connected directly or via a V-belt, (this system uses a V-belt: its 
design will be discussed shortly). The generator is designed with a built in voltage regulator, 
however in system testing, it was found early on that that regulator was not robust enough to 
control the voltage through the changes in load that the system required." 1 

Throughout this project, the exact 
operating characteristics of the generator 
have been unverified. The primary reason 
for the choice of generator was price. The 
Ming Dong generator was manufactured in 
China and has very few published 
specifications; those that are published can 
be found in Appendix A. The unknown 


Figure 3: Adjustable test-bed of loads: This 
arrangement of loads from 
75W per phase to 1125W was created for initial 
testing. 
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nature of the generator makes the effectiveness and reliability of the feedback system, as well as 
system testing, incredibly important. 

The load used in this model is a Y-connected arrangement of light-bulbs. Only these 
resistive loads have been used to test the model system so far. Each line on the Y-connected 
load has three inserts for light-bulbs that can be switched on or shorted; this arrangement 
provides a discrete variable load for the preliminary testing done in this project (see Figure 3 for 
image of the load arrangement). A complete parts list can be found in Appendix B. 

3.2 Sealed Control System 

In this system there are two independent control loops that mirror the governor and 
voltage regulator of the DDG plant. They control the frequency and voltage RMS values, 
respectively. 

3.2.1 Governor 

DC motor speed control can be implemented in any of three ways, via field resistance 
control, armature voltage control, and armature resistance control. xvn This design controls the 
motor via armature voltage with proportional and integral control (PI control). With just 
proportional control, the system feedback is more likely to overshoot the desired frequency 
beyond the acceptable range of the Navy specifications. Therefore, integral control is also 
needed. However, the less the relative weight of the integral control constant relative to the 
proportional control constant, the faster the system will return to 60 Hz after a major 
disturbance. xvl " 

Ideally, the speed of a DC motor is directly proportional to the input DC terminal voltage. 
In reality, resistance to increases in torque for upon increased terminal voltage is caused by the 
back emf, this reduces the slope of the voltage-speed curve at higher voltages. xlx The most 
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common means of controlling a DC motor is to use pulse-width modulation control to vary the 
terminal voltage." However, for the sake of expediency and simplicity, this system uses a 
remotely controlled power supply to generate the terminal voltage to the motor at a level 
calculated by the software based feedback system. 

The feedback loop measures the speed of the DC motor against a reference frequency, by 
default 60Hz. While the actual frequency of the system is from the output of the generator, that 
value is mechanically tied to the speed of the motor and the turns ratio of the V-belt. A Hall 
Effect sensor on the output of the motor is a straightforward and acceptably accurate means of 
reading the speed of the motor, and with equation 1, the frequency of the system. 

The software control loop sets the power supply remote control voltage, V cr , between 0- 
10 V. This value controls the output of power supply, ie, V tenn of the DC motor. Therefore V cr 
directly controls output speed and V cr adjustments can remove the error between the DC motor 
speed and 60Hz. A diagram that includes this loop is shown in Figure 4. 

3.2.2 Voltage Regulator 

The second control loop of the scaled 

system provides feedback on RMS values of 
the output phase voltage of the generator. The 
voltage regulator for this system is 
implemented as shown in Figure 4. A PI 
control loop was chosen for the feedback 
system. The voltage regulator for this system 
is based on the error between the sum of the 
three phase-to-ground RMS voltages and a 


Figure 4: Control Loop Diagram 
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reference voltage. The feedback loop directly adjusts the field current of the generator rotor. 

The field current is adjusted with a variable power supply that sends a DC voltage through the 
field resistance (again, a power supply was chosen expediency and simplicity). 

The actual control loops are implemented in the Lab VIEW development environment. PI 
control loops were chosen as the control method for both the governor and voltage regulator 
loops given the popularity of that method of control for motor-generator control loops as 
published by Stephen Limans’™ and Kim et al xx “. 

The details of the signal manipulation will be discussed in the next section. The design 
of this system was not chosen for power consumption efficiency or practicality in applications. 
The purpose of the design is a proof-of-concept scaled power generation system for a future 
laboratory test-bed for the load monitoring system. Therefore the choice of parts for the system 
was driven by the team's best judgment of what would be effective but the fastest to build and 
cheapest. 

Scaled System Operating Requirements 

The performance requirements for the scaled system will be the same as those specified 
for the destroyer system. Given the responsiveness of the DC motor system compared to a gas 
turbine system, one concern was that the scaled system would be too responsive to changes in 
load to be an accurate model for the gas turbine generator. In this case, a time delay would have 
to be added to slow the feedback system down. Tests, as described in section 6, were done to 
measure the responsiveness of the system. The physical differences between the scaled system 
and the naval system are summarized in Table 1. 
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Table 1: Comparison of Requirements of Naval and scaled systems 



Scaled Model 

DDG 

Prime Mover 

Two parallel DC motors 

Allison 501-K34 Gas Turbine Engine 

Generator 

STC-5 series 3-phase synchronous 
Mingdong 

Allison 501-K34 Gas Turbine Generator set 

Rated Voltage 

400V/230V 

450 

Rated Current 
(A) 

9 

4009 

Power Factor 

0.8 (lag) 

0.8 (lag) 

Output KW 

5 

2500 

Rated Speed 

1500RPM 

14.340 RPM 

Pole Number 

4 

2 

Exciter 

DC power supply 

Three phase brushless excited motor 


4. Design of Scaled Model 

4.1 Software Modeling 

In order to prove that the DC motor-generator would suffice as a model for the gas 
turbine-generator system it was necessary to model the two systems in a software modeling 
environment. The modeling work and simulations were completed in the theses of two U.S. 
Navy engineering duty officers, Gregory Elkins xxm and Jeremy Leghorn.™' They developed 
models of the motor operation, LCDR Elkins, and the generator, LCDR Leghorn, to prove the 
system could be used as a scaled DDG model. 

A per-unit system was developed in MATLAB Simulink that models the behavior of the 
physical scaled system in the lab. Power per unit is 5KW as that is the rated real power for the 
generator. Due to the ratio of the V-belt between the DC motor and generator, the speed per unit 
is 2880RPM. This motor speed runs the generator at 60Hz. 
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4.2 Detailed Hardware and Software Implementation 

The following section discusses in detail the design choices of the current physical 
feedback loop. Work on the implementation of this system was done through UROP 
(Undergraduate Research Opportunity Program) by Jacob Osterberg, and by the author for this 
6.UAP (Undergraduate Advanced Project). Circuit diagrams and further details can be found in 
Appendix C. 

The design priorities on the implementation were cost and expediency. Therefore when 
possible, equipment available in the laboratory and what we thought would be quick, reliable 
solutions were chosen. In some cases, the actual expediency of the design choices was less than 
expected. This section includes a discussion of lessons learned and implementation challenges 
as well as the final design. 

4.2.1 Software and I/O Card 

The first purchase and defining element of the feedback control project was the 
software/hardware interface for the control loop. The interface built consists of a hobbyist kit 
input-output board and the LabVIEW Student Edition software development environment. A 
software-based control loop was chosen over pulse-width modulated power electronics with the 
belief that it would provide more flexibility and extensibility for future testing. Hopefully, this 
system will be the foundation for the streamlined user interface of the future TableTop test-bed. 
The multifunction input-output card is the connection between system control signals and 
measurements, it is controlled by the LabVIEW user interface. The card that was chosen has a 
total of thirty-three digital and analog inputs and outputs, although this is more than necessary, 
this board is large and flexible enough to incorporate future additions or testing regimes without 
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additional investment. Additionally, with all inputs and outputs controlled from one central 
location, data acquisition and analysis will hopefully be easier. 

The user interface for the physical model was built in LabVIEW. The primary LabVIEW 
Virtual Interface (VI) is called the "Panel.” The Panel controls and reads physical values 
through the input-output card which is a Velleman K8061 USB Extended USB interface card 
(also known as the preassembled VM140 card). The card will be referred to as the I/O card 
(shown in Figure 5). The TO card comes with a software driver and a dynamic link library 
(DLL) with functions that command the TO card. In LabVIEW, a LabVIEW Library was 
created as the central location for all functions and code relating to this project, the library is 
called the User Interface Library (UILib). The I/O card DLL was imported to UILib and every 
card function (there are twenty-five) was compiled as an independent sub-virtual interface 
(subVI) that can be run individually within other Vis in UILib. The user manual for the K8061 
card can be found online or on the laboratory computer. The Panel was created as a VI in UILib 
as well; the function subVIs are called from the Panel. 

Figure 5: I/O Card. The input/output card for the system 
was a K-8061 Velleman Hoobyist Kit USB Board. 



Although the LabVIEW-I/O card arrangement provides flexibility and user-friendliness, 
there are limitations to the I/O card’s capabilities that required an unforeseen increase in the 
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complexity of signal manipulation between the system and the I/O card. The I/O card has both 
analog and digital inputs and outputs, however, it translates all signals to either 8 (for output) or 
10 (for input) bit digital values on either a 0-5V or 0-10V scale. Precision is lost. Furthermore, 
the channels are specified as DC voltage channels by the manufacturer. The reason for this 
specification is that calls to read or write to a channel on the card are discrete and only return or 
set a voltage one time when they are called. These specifications for the I/O card have required 
increased complexity in software and hardware for this loop. 

The hardware adjustments necessary for use of the I/O card were to translate both the 
motor speed (approx. 2880 RPM) and the phase-to-ground RMS voltages (approx. 120 V) to DC 
voltages between 0-10V that can be read into the I/O card. The circuits for both of these 
translations were built primarily by Jake Osterberg (see Appendix D). An 8-bit PIC18F4620 
microcontroller was used to translate the motor speed to a DC value and LTC1966 Precision 
Micropower RMS-to-DC Converter chips were used to read the RMS voltage levels of the output 
phase voltages. The chip translates the RMS voltage to a DC value between 0-10V. Although 
the addition of these circuits has added complexity to the system, complexity itself is not a 
problem. Flowever, it has increased the likelihood that a small error or noise in any part of the 
system will magnify and propagate to the voltages read by the software. To counter the 
uncertainty added by these parts, control system constants have been adjusted. The frequency 
reading is very consistent during steady state. Figure 8 shows stead state frequency readings. 

The discrete jumps are the difference between the 8-bit digital values that voltage is converted 
into from the I/O board. 
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Figure 6: Steady state frequency during feedback. The variations are discrete and 
caused by rounding between 8-bit input values from the I/O board. 



4.2.2 Power Supplies 

Sources of DC voltage for this system are needed to power the DC motor and the field 
current in the generator. We decided to use power supplies for both of these sources as they 
would be the simplest to implement from the board as long as the supplies have analog voltage- 
based remote control and provide enough power. The initial goal for the optimal terminal 
voltage power supply for the motor voltage was 200V and 10A. However, when two 150V, 7A 
supplies that could be paralleled to provide 150V and 14A were found in the lab, they were used 
instead. xxv Given that researching the price for a controllable supply with these high power 
ratings made it clear that one would not be purchased for less than $2000, these power supplies 
were the best choice. 

These supplies are 1000W series XHR 150-7s made by Xantrex (see Appendix E). The 
power supplies are remote controllable and when placed in a master-slave configuration, they 
should produce 14A and 150V. This power output fits the specifications for system power with 
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an adjustment in the motor to generator pulley. XXVI The lower power at the input of the system 
requires a lower terminal voltage to be able to drive the generator at rated frequency, 60 Hz. 

The pulley adjustment allowed a lower motor speed to correspond with a higher generator speed. 

The XHR 150-7 output can be remotely controlled with a 0-5V or 0-10V DC voltage that 
will linearly correspond to an output voltage of 0-150V. We set the system to run on a 0-10V 
DC input from the I/O card. Each discrete 
voltage level specified by the 8-bit I/O 
card analog output is a jump of .58 volts 
on a scale of 0-150V. The implications 
of this constraint are discussed in the 
following section. The voltage-speed 
curve of the motor is shown in Figure 6. 

The field current power supply 

was purchased very late due to the high cost of the supply needed. It is calculated that we will 
need a power supply that provides simultaneously 3A and 50V. Currently, there is manually 
controlled power supply that fits these specifications in use in the system. The purchase delays 
were because no analog signal control supply could be found for less than $1000. However, an 
economical supply operated by serial port was found, and has been purchased though not 
incoiporated into the system. It is an Instek PSP-603 60V/3.5A supply. This supply can be 
integrated into the LabVIEW Panel using the control functions through the USB serial port 
instead of analog control through the I/O card. 


Figure 7: Speed-Voltage Curve ot Motor 


Motor Constant 
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4.2.3 Pulley System 

The physical connection between the motor and generator is a V-belt. The sizing of the 
pulley between the two is essential to meeting the operating requirements of the system under 
loads given the restrictions of power from the XHR 150-7s. Initially, the pulley had 32 teeth. To 
calculate the speed requirements for this system we can use equation [2]. We know our 
maximum voltage and current is 150V and 14A respectively. Also, from characterizing the 
motor we know K and R from the voltage-speed plot shown in Figure 6. The pulley was 
adjusted to one that has 22 teeth so that rated speed would be reached at approximately 125 V 
unloaded. 

4.2.4 Summary of System Operation 

The first step in the system is when LabVIEW (either through automated feedback or 
manual user control) designates a terminal voltage and, in future, a field current level that will 
drive the motor and generator, respectively. These voltage levels will be translated to 8-bit 
values and then outputted as the corresponding voltages between 0-10V on the I/O card. The 
terminal voltage 0-10V DC voltage is inputted to the back of the master XHR 150-7 power 
supply and commands the output voltage of the power supplies. 

The terminal voltage is applied to the armature of the DC motors and produces 
mechanical torque in the series motors. The shaft of the DC motors is connected via V-belt to 
the rotor of the generator. Mechanical power is transferred to the generator via V-belt. The 
speed of the motors determines generator frequency. In a separate loop, the field current 
determines generator output voltage amplitude. In order to control the terminal voltage and field 
current, the frequency and RMS system voltage are feedback to the software to calculate the 
appropriate next value for these signals. 
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5. System testing and Analysis 

The graduate theses on this system were to test whether the physical system, as modeled 
in software simulations could theoretically represent a scaled version of a DDG gas turbine 
generator. The purpose of this undergraduate project was to prove or disprove the physical 
ability of the system to operate as a scaled version of the DDG system. The answer to whether 
this physical system can be used in future testing as a reduced scale version of the DDG plant is 
yes, with a caveat. The present arrangement has major limitations in two areas, some can be 
improved upon, some can not. The first source of limitations is the power supply rating. The 
second problem that hinders the successful replication of the successfully simulated model in the 
physical system is the imprecise nature of both the input and output signals of the control loop. 

In testing, the control loop constants had to be adjusted to fit the time and accuracy limits 
of the physical system. The optimal control constants for the physical system are very different 
than those of the Simulink system. One major difference between the two is that the physical 
system has a delay time between iterations of the feedback loop. This component was not 
considered in the more complex mathematical models in Simulink. The equations used for 
feedback in the physical system are shown below. The variables are defined as follows, e is the 
error, fset is the desired system frequency, f sy stem is the frequency read in for the system. 


fsit fsystem 


(4) 



Although it is not in use now, the secondary loop, for RMS voltage and field current is 
set up in LabVIEW. If the inheritors of the project have access to a 0-10V remote controlled 
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power supply to power field current, all that will be necessary is a tuning of the gains in 
LabVIEW. A serial-port controlled power supply will require additional work simply to add the 
control methods to the UILib and call those instead of I/O card DLL functions to set the voltage. 
Data acquisition is already running for RMS voltage and will be able to run through the board as 
set up. Another channel can be used to record the field current power supply voltage if desired. 

In development of this project, major errors were recorded during testing related to 
magnetic field interference. These errors produced erratic frequency readings. An example of 
the incorrect data is shown in figure 7. After debugging, the problem was discovered and 
physical layout changed. The changes removed the error. After that improvement, tuning could 
take place. 

Instability invariably 
occurred when the physical 
system was put under the control 
conditions found in the Simulink 
model. The reason for this 
instability is that the simulated 
governor is built to correct 
frequency to within 0.15 Hz. In 
the physical system, the use of 

the relatively large PI constants developed for this precise control causes extreme oscillation. 
Using these constants demands the system operate more precisely than it can measure. The 
feedback control is trying to maintain an error less that 0.15 Hz and the speed reading is 
oscillating occasionally more than that. See figure 8 for an example of the steady state error in 


Figure 8: Data Run with Data Errors. Red is where unexpected and 
extreme frequency values were read into the system. 



97 














20 

frequency. The error is mainly due to rounding errors of the 8-bit digital representation of the 
frequency in the software. This limitation is caused by the I/O card. Instead of stabilizing 
frequency, with this error, the simulated constants cause the governor to alternately force the 
terminal voltage to zero and maximum, 150 V. It was found that an initial K p = 0.01, K; = 0.0, 
and feedback loop period = . 1 sec could control the system stably and near Navy regulations. 
These values are very far from the Simulink values of K p = 1.8, K; = .23 and feedback period = 
nearly continuous. 

After tuning, the optimal constants were found to be K p = .0012, Kj = .0012, time delay = 
.08 seconds and no dead band. The results from a load change from 225W to 675W is shown in 
Figure 9. The relationship between the K-values and time delays is not mathematically known 
and needed to be tuned manually. It was roughly found that the faster the feedback operates, the 
lower the K-values needed to be to avoid overshoot. The optimal response is the best fit for the 

Figure 9: Good Feedback Data. Example of good feedback occurring under a change of load from 

75W to 225W per phase. 
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Simulink model because it doesn’t overshoot and it also smoothly returns to steady frequency 
like the GTG system. 

A major limitation of the current system is that the power supply simply can't deliver 
high enough power for full-load operation. When the maximum load of the current load 
platform. 1 125W, is applied the system saturated. The effect on the frequency can be seen in 
Figure 10. This figure shows an application of the 1125W load at approximately 50 iterations 
and 150 iterations. The response to the application of the load, as well as observation of power 
supply operation, show that the power supplies were current limited during the response. In 
order to handle rated load for the system, actually 5kW, a much more powerful power supply is 
needed. The actual response of the GTG system is shown in Figure 1 1 for comparison. 101 ™ 


Figure 10: Largest Applied Load. 1125W applied and removed from system with Ks 

= .0012 and time delay = .08 sec. 

Application and Removal of 1125W 



Another place where imprecision is added to the system is where the I/O card sets 
commands the terminal voltage level. Output analog voltages from the card are 8-bit values. 
Every step in the 8-bit value corresponds to a change of .58 V in the terminal voltage. This 
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Figure 11: Response of the actual GTG system to application and removal of its rated 



problem, though not conveniently solved due to the limitations of the board, can be improved. A simple 
improvement would be to split the 0-10 V power supply remote control signal between two analog 
outputs between 0-5 V and then use an operational amplifier to sum these two voltages. The resulting 
terminal voltage control signal would be precise to 512 partitions of 0-10 V instead of only 256. This 
circuit is built and is an easy extension for the project if more precision is found to be necessary. 

Despite the improvement that can and should be made to this system, it is currently set to operate 
as close as possible to the GTG operating characteristic. For instance, when a smaller load that does not 
saturate the power supplies is applied, the model's response is the closest it can be to the GTG model. 
The response shown in Figure 12 is similar to the actual GTG model. The response of the scaled model 
was tuned to fit this response as best as possible. 
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Figure 12: The response of the system to less than rated load (225W to 675W and back to 
225W). This is the closest the scaled system can come to the GTG model. 



6. User Interface User Manual 

The user interface for the model is implemented in LabVIEW. Front panel testing can be done 
easily for manual, step, or variable feedback control. The Panel is shown in Figure 13. Operation from 
this screen is straightforward. However, editing the code for the interface is much more complex. Both 
operations will be discussed in this section. 

6.1 Panel Control 

The Panel is designed to provide some flexible testing and control options. Iterations of 
testing for the software and for the system models drove an interface design that has multiple 
operating methods, the user must designate the operating method prior to running the control 
software. There are three options that are available for operation; 
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Figure 13: The Panel User Interface 


MANUAL CONTROL without FEEDBACK 


INSTRUCTIONS: 

Before pressing play, designate "Type of Run” 
and "Total Time of Run (sec)." (Save data input 
is not necessary, it will be prompted at the end of 
the run.) 

If you are running feedback, enter K-values pnor 
to starting run. 

Type of Run 

^[Manual Control (no feedback) 

Run (sec) Card Addr 


Armature Voltage (0-150V) 



Vcontrol (0-10V) 

Jo :jo.o5 

Digital Input (0-255) 

0 


Delay Time Manual 


| Save Data to File in Manual Control 
(promt after run) 




OPEN LOOP STEP RESPONSE 


Step (sec) 

IT- 


Time elapsed since step 
(T“ | STOP I 

Start Voltage (0-150V) Voltage After Step V dig 

o : : Jo Jo 


Elapsed Time (s) Card connected? 


MANUAL START, transition to FEEDBACK (feedback runs for "run" time input) 


Vcontrol (0-10V) 


Save Data to File 
(promt after run) 


Digital Input (0-255) 



Armature Voltage (0-150V) 
75 

50 i 

25- J “125 


100 



Delay Time (sec) (~T2) 



Field Current Control: Equation is: (Kp*error+integral(Ki*error)] * Vail + 
Will control around sum of three phases 




1. Manual control of DC motor input voltage and generator field current without 
feedback. This set-up was designed for preliminary testing. It remains on the 
user interface so that it is easily accessible if new elements are integrated into the 
system and manual control is necessary for testing. 

2. Open-loop step response of motor input voltage provides an easy set up for open- 
loop step response testing. 

3. A feedback system. In this option manual control of V cr is allowed for an 
undesignated amount of time so that the operator can put the motor in whatever 
steady state condition they wish to start the test in, operations during this time are 
not saved or timed. When ready, the operator presses the ‘Start Run' button and 
the feedback loop takes over V cr control. The feedback loop will run for either 
the pre-designated run time or until the ‘Emergency Stop’ button is pressed. 
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6.1.2 Running 

After opening Panel and before hitting the LabVIEW run button, the user must designate the run 
type. They also must designate the total run time. For open-loop step response and control, further test 
variables should be set before the program is run. However, these values are adjustable during operation 
- especially for feedback control prior to pressing the ‘Start Feedback' button. 

When the runs LabVIEW, the Card Address should change to 0 and ‘Card Connected?’ should be 
one, representing true. If this is not the case, the USB is not plugged in, or there is a major problem with 
the block diagram. The next step will likely to be to adjust the voltage manually (for manual and 
feedback control). Voltage input can be controlled by dial or by typing into the textual prompt to the left 
of the “Armature Voltage’' dial. The input value to this controller will be exactly the output of the power 
supply. 

For feedback, there are important settings to pay attention to. The K p and K, values as well as the 
time delay need to be set on the Panel. Also, the dead band size will by default be .08 multiplied by the 
average of the gains and then by 150. This can be adjusted in the block diagram if desired. 

6.1.3 Saving 

Data is recorded at discrete time intervals from every analog input on the I/O board. The 
saving time interval is the same as the feedback interval during feedback, it can be set on the 
Panel in the other settings. At the end of a run, the Panel will automatically prompt the user for a 
file path to save the data. The data is saved in vertical columns. Table 2 provides a list of the 
column headers. The Panel saves all analog inputs to the files it creates. Even though some 
columns are empty under the current set-up, this saving method provides allows for easy 
extension of data acquisition. The data can be imported to excel or MATLAB for analysis. 
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Table 2: Column Headers in a Panel data acquisition file 


COLUMN 

MEASUREMENT 

RANGE 

/UNITS 

1 

Time 

Time Delay 
(seconds) 

2 

System Frequency 

50-70 Hz 

3 

Terminal Voltage 

0-150V 

4 

Vrms 0)a 

0-250 V 

5 

Vrms 

0-250 V 

6 

Vrms 

0-250 V 

7 

Frequency Error 

Voltage 

8 

Integrator Sum 

Voltage 

9 

Dead Band Indicator 

Boolean 


6.2 The Lab VIEW Block Diagram 

The code for the Panel is graphically designed a block diagram. A snapshoot of part 
of the block diagram is shown in Figure 14. Future operators of the system can adjust the design 
of the control system by editing this code. The basic format is a sequential loop that commands 
the I/O card. The card must first be opened before it can be controlled. Then, depending on the 
user’s choice of operating scheme, an operating loop, usually a while loop, will execute until the 
time is expired or a stop button is pressed. As the while loop executes, board data is read into an 
array that builds indexes every time the loop executes. At the end of the loop, the data array is 
saved. After saving, the card is cleared of data commands and then closed. Sub-Vis that execute 
data conversions, control methods, and the actual feedback have been designed to simplify the 
highest level loop for the system. 
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Figure 14: Part of the LabVIEW Panel Block Diagram. LabVIEW Panel has many 
layers of block diagrams to simply the graphical code. 



The I/O card commands are in the form of sub-VIs and are stored in the ‘Vis’ folder 
within the UILib. Depending on whether the operating system of the computer uses 32 or 64 
bits, there are different DLL libraries that must be used. To switch between the two, one must 
recompile the library and replace each sub-VI. There is also a small but major bug in the DLL 
library in the 'ReadAnalogChanneL function. This essential function does not return the data in 
as written in the DLL library, however, the operator can fix this problem by editing the header 
file and ensuring that the data is returned manually during the LabVIEW compilation of the 
functions. 

7. Conclusion 

The results from this project show that a future Table Top test bed is most likely 
possible with the addition of more powerful power supplies. It may be possible to more closely 
approximate the GTG response on this system if these supplies were larger. Also, improvements 
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could be made in the control loop with a possible inclusion of a derivative term for a smoother 
return to steady state. Maximum load could not be applied in these samples due to both the 
limitations of the power supply and the load test-bed. However, Figure 13 shows a summarizing 
image of load application and removal with working feedback. The load changes from 225W to 
675W in this run. 

Many lessons were learned in this project. The use of LabVIEW was central to this 
project and the programs strengths and weaknesses became clear to the author, strengths being 
power and depth of the National Instruments systems, and weaknesses dominated by the 
inconvenience of not using National Instruments hardware and Toolboxes for implementing the 
interface. Despite this, I still believe that LabVIEW has provided a very user friendly user- 
interface and easy incorporation of the third party DLL library. 

Although magnetic interferences delayed the project and testing for many days, the 
effect of this phenomenon was a good lesson to have hammered home. The first few reactions 
were that the problem was in hardware, software, or the generator itself. However, it is 
important to always check physical connections and influences on the system. These major, 
game-stopping problems were quickly fixed once the cause was found. 

In order for this problem to move foiward, the software will be placed on a computer 
in lab. The drivers for the I/O card and LabVIEW will be stored in the lab. Also, for the aid of 
those who continue this project, the appendices of this paper will include more details on 
hardware and software operation. A file will be saved to the lab computer that includes these 
files as well as additional background information. 
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Appendix C: Additional Generator Information, by Vanessa Esch 

Appendix A: Additional Generator Information 


Generator Specifications were found at the following website from the manufacturers in China. 
The generator was bought off eBay. 

i mi>iny Contact 

STC generators arc to be used in town, the countryside, worksites, mountain and pasture lands as a elertric 
power source Tor sighting purpose, It car also be used as a reserved power source Tor emergert case. 

The generators a^e of d ri p-poof with rotary field type and with tile adoption of harmonic excitation 
system, that alow your easy operation and simple maintenance. The generates a*e of three-phase four-wire 
type, making use c* alar correction with neutral point. The rated Sne voltage s 400V. phase voltage 230V. 
frequency 50Hz. power factor 0.8{lag). It can provide 60Hz and the other voltage s generator according to 
necessary. 

They car be coupled with a prime mover drectfy or through V-helt making right or reverse continuous 
rotation at the rated sceed. when the revolution slip of prime mover is 3% and load varies in the range of 0 
100% cos 1.0-0.8 they give feed constant vc tage. After sudden charge { -crease or decrease’lof load, the 
generators w9 scon return to their normal working state, without ary starting devices the generator car directly 
start an unloaded squirrel cage induction 
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Appendix D: Parts List, by Vanessa Esch 


Appendix B: Parts List 

Terminal Voltage Power Supply: (2 x) Xantrex XHR 150-7 Power Supply 

Field Current Power Supply: Instek PSP-603 Programmable DC Power Supply, 60V/3.5A 

Input/Output Card: Velleman K8061 Hobbyist Board (also called, preassembled, a VM140) 

DC Motor: (2 x) Leeson Direct Current Permanent Magnet Motor. Model: C42D340T7A, 
Volts: 120V. RPM: 4800 RPM, Current: 14A, 2 H.P. 

Pulley: McMaster-Carr Drive Pully, L-series, 22 teeth, Vi\ Bushing bore. Part Number 
6495K213 

Generator: STC-5 Ming Dong Generator, 5KW, 

Speed-to-DC Converter: (See Appendix C) 

RMS-to-DC Converter: (See Appendix C) 


4 pole machine 

Note: Converters are shown in 
Appendix F of this thesis. Addition 
to the parts list made below. 

- Greg Elkins 


Addition 


Parts listed below were purchased through McMaster-Carr (www.mcmaster.com) unless 
otherwise specified. 

5KW 3-Phase Generator: listed above, Purchased through Georgia Generator on ebay 
(stores.shop.ebay.com/Georgia-Generator) operated by Tom Osborne (478-457-5524) 

Motor Pulley: listed above 

Motor Pulley Quick-Disconnect Bushing: Style JA, 5/8" Bore, Part # 6086K112 

Generator Pulley: Timing Belt Pulley, L Series, 32 Teeth, OD 4.06", Bushing Bore, 

Part # 6495K218 
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Generator Pulley Quick-Disconnect Bushing: Style SDS, 1-1/2" Bore, Part # 6086K324 



Timing Belt: Trapezoidal Tooth, Neoprene, L-Series, 1/2" Width, 39" Outer Circle, 
3/8" Pitch, Part # 6484K157 

Motor Coupling (connects 2 motors): Steel One-Piece Set-Screw Coupling, 5/8" Bore, 
2" Length, 1-1/4" OD, Without Keyway, Part # 6412K15 



Gearbox (for future consideration): Howse Gearbox for M60 Rough Cut Mower, 55HP, 
1.47 Gear Ratio, Model # 45PRC30147-N, available at Northern, Tool + Equipment 
(www.northemtool.com) 
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Appendix E: Pictorial, by Vanessa Esch 


Appendix C: Pictorial Walk-Through System 

1. The first step in the control of the system is setting the Terminal Voltage and Field Current from 
the Lab VIEW software interface. The laptop is shown as arranged during testing. 



2. A USB connects the laptop and K8061 card. 



3. The output from the K8061 card and ground go directly to the back of the Power Supply (red and 
black) to command its output voltage. Current out of the terminal voltage power supply is not 
actively controlled but is a function of the load on the system. 


Ill 


















4. Terminal voltage from the power supply is applied to the prime mover DC motor rotors in series. 
(Yellow wires) 



5. A Hall Effect Sensor reads the speed of the prime movers. (See previous picture, far side of 
motors.) 


6 . 


The digital output of the Hall Effect Sensor is sent to a conversion board that outputs an analog 
voltage between 0-10 V that corresponds to the frequency of the system. Board on the left. 
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7. The prime mover mechanical torque is transferred to the generator via V-belt. 



8. The output voltage from the generator goes to the load bed 
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9. The output from the generator is stepped down and then applied to the RMS-to-Analog 

conversion board. The output of this board is set to the K8061 board and then to the software 
feedback loop. 
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Appendix F: Signal Conditioning, by Jacob Osterberg 


Generator Feedback Electronic Hardware Specifications 

Jacob Osterberg 
Spring 2009 

Overview: 

Two additional circuits were necessary beyond the Velleman Extended USB 
Interface Card in order to interface the generator system and the measurement computer. 
Both speed and voltage measurements of the system were necessary. Since the extended 
USB interface card did not have input hardware interrupts, it would have been difficult to 
accurately measure the TTL input pulse from the Hall Effect speed sensor on the motor. 
Therefore the frequency converter circuit was created to convert this TTL pulse into an 
analog voltage easily readable by the USB interface card. The generator voltage 
measurement may have been possible through fast sampling with the USB interface card, 
but an outside RMS conversion would be less taxing on the card, provide a more accurate 
measurement, and reduce the software complexity in Labview. For this reason a RMS 
converter circuit was built to measure, amplify and filter the generator voltage. Finally a 
modified version of the frequency converter was designed as a PCB with some additional 
capabilities. 

RMS Converter: 

Function: 

This circuit receives an AC voltage waveform with an amplitude of 
up to 350V and converts it into a RMS value scaled from 0 to 10V. 

Details: 

Three identical circuits (one for each phase voltage) use a resistor 
divider to scale down an AC input waveform from a maximum amplitude 
of 350V to a IV maximum AC input to the Linear Technology LTC1966 
RMS chip. Since the LTC1966 does not have a negative voltage supply it 
was necessary to insure the input AC waveform does not fall below 0V. 
This is accomplished by using a simple resistor divider creating a 
reference voltage of 2.5V for the IV AC waveform to oscillate around. 

The output of the LTC1966 is then averaged by a luF capacitor on its 
output. Next a passive filter with a cutoff frequency of 15.4 Hz further 
smoothes the waveform. Finally the signal is amplified by an AD358 op- 
amp resulting in a final output voltage of 10V for a 350V peak AC input. 
As currently configured a RMS voltage of 208V corresponds to an 8.2V 
output from the circuit and a RMS voltage of 120V corresponds to a 4.8V 
output. 

Discussion: 

When selecting resistors for the circuit it is essential that great care 
be taken in measuring and noting the actual value of the resistances and 
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not their stated value. Most resistors have tolerances of 5-10%. This 
variation of 5-10% can greatly affect the resistor dividers and the op-amp 
amplification circuit. While it is possible to purchase resistors with 
smaller tolerances, the higher tolerance resistors can work fine. Through 
direct measurement, only resistors with values close to the desired value 
will be selected. Also in the final circuit calculation actual measured 
values should be used instead of the ideal printed values on the resistor. 

An averaging capacitor of luF was selected for the output of the 
LTC1966 RMS chip as recommended by the datasheet. An additional 
passive filter with a cut off frequency of 15.4 Hz was also added to 
smooth the output waveform. These act as a significant delay to the 
voltage measurement, but greatly reduce the noise. The passive filter may 
have to be changed for future measurements depending on the acceptable 
delay and noise when using the RMS signal. 



Figure 1: RMS Converter Schematic 
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Figure 2: RMS Converter Circuit Photo 


Frequency Converter: 

Function: 

This circuit receives a TTL square wave pulse input and converts 
the frequency of that pulse within a specified range to a 0 to 1OV analog 
output with 256 discrete values. 

Details: 

A TTL square wave input ranging from 40Hz to 2500Hz is input 
into a 74HC14. The signal is then sent to a PIC18F4620 microcontroller. 
Depending on constants set in internal code the microcontroller outputs a 
scaled 8 bit value to an AD558 DAC. The result is an analog voltage 
which corresponds to the input frequency. The maximum output voltage 
is 10V with a minimum of 0V. There are 256 possible steps in the analog 
output voltage. 

Discussion: 

Currently the input pulse is received from a Hall Effect sensor 
connected to the motor driving the generator. The sensor triggers 14 times 
for every revolution of the motor. At an output frequency of 60Hz from 
the generator this corresponds to a 611Hz signal from the sensor. This 
circuit is set to output 0V from 0-509Hz (there is a slight glitch at 
38.147Hz from the internal timer overflowing, but it has a very small 
<0.1Hz bandwidth and should not effect the general operation of the 
circuit). The sensor frequency of 509Hz corresponds to a generator 
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voltage frequency of 50Hz. At 611Hz (a generator voltage frequency of 
60Hz) the output will be 5V and at 713Hz (a generator voltage frequency 
of 70Hz) the maximum output voltage of 10V will be reached. All 
measurements from the sensor at greater than 713Hz will result in a 10V 
output. 

The frequency conversion can be easily changed in the PIC 
software by modifying the constant values in “frequency.h”. This allows 
for an easy change from 509-713Hz to anything desired within the 40- 
2500Hz range of the current setup. Since there are only 256 values 
available from the DAC this range should be modified to suit the needs of 
the system. A larger lObit DAC could also be substituted with minor 
additional programming in order to increase the resolution of the system to 
1024 values. 

The 74HC14 Schmitt Trigger may seem redundant in this circuit, 
but it provides a buffer between any inputs and the microcontroller. Also 
if the input signal is not perfect the 74HC14 will ensure a nice edge 
delivered to the microcontroller interrupt. The 74HC14 also has an 
oscillator circuit built on it with a potentiometer in order for quick circuit 
testing to be preformed without the aid of a function generator. 

The PIC18F4620 is overkill for this application, but it was familiar 
and available at the start of the project and provided the easiest 
microcontroller option. It also gives future flexibility for the entire 
feedback control loop to be moved into the microcontroller. For 
programming, I recommend using the PlC’kit 2 

(http://www.microchip.com stellentidcplg°IdcService=SS GET PAGE& 

nodeId=1406&dDocName=en023805) . It is a cheap programming option 
for PIC microcontrollers. For software, I used the MPLAB ICD software, 
available free from Microchip 

0ittp://www.microchip.com stellent idcplg?IdcService=SS GET PAGE& 

node!d= 1406&dDocName=en019469&part=SW007002 ), with the free 
student MPLAB C18 C complier also available from Microchip 
(http://www.microchip.com stellent idcpl 2 : 0 IdcService=SS GET PAGE& 

nodeId=1406&dDocName=en010014&oart=SW006011 ). 
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Figure 3: Frequency Converter Schematic 



Figure 4: Frequency Converter Circuit Photo 
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Frequency Converter PCB: 

Function: 

Tins printed circuit board functions very similar to the previous 
frequency converter built on the prototyping breadboard, but it supports 
additional inputs and outputs including a serial port. 

Details: 

The same basic hardware is all present from the original frequency 
converter circuit. It has the ability to sense a TTL input and output a 
corresponding analog voltage. Additional features include access to more 
inputs and outputs and the ability to communicate to a computer through a 
serial port. Also switches and LEDs have been added for user 
functionality. Although all these features may be useful, it is not 
necessary to populate the entire board for it to perform the frequency 
conversion function. The MAX3221. DB9 and supporting capacitors are 
not required for the frequency function of the circuit. Take care to ensure 
that the correct pins are connected (by a jumper connector) on the AD558 
for either 2.56V operation or 10V operation. When programming ensure 
that the white triangle on the PICkit 2 programming port lines up with pin 
1 on the connector (also pin 1 on the PIC18F4620). 

Discussion: 

The functionality of this version of the frequency conversion 
circuit is much more than a simple measurement, calculation, and output. 
This circuit is capable of containing the entire feedback loop for the 
generator system, hi addition, data can be sent to and from a computer 
connected through the serial port. The RMS voltage may also be 
measured by a built in lObit analog to digital converter in the 
microcontroller. The possibilities of this circuit are also not confined to 
the generator system. This board provides many inputs and outputs useful 
for additional applications. When used with the PICkit 2 programmer 
from Microchip, this board becomes a very powerful tool to quickly 
measure and respond to many different signals. 
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Figure 6: Frequency Converter PC'B Board 
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Parts List for Frequency Converter PCB: 

Cl 1 uf Through Hole 

C2 47nF Through Hole 

C3 .01 uF Through Hole 

C4 .OluF Thiough Hole 

C5 luf Through Hole 

C6 .luF 603 Surface Mount 

C7 . luF 603 Surface Mount 

C8 . luF 603 Surface Mount 

C9 .luF 603 Surface Mount 

CIO . 1 uF Thiough Hole 

D1 Red LED 

D2 Red LED 

J1 2.1 min Jack 

J2 5 Pin .1" Spacing Header 

J3 6 Pin . 1" Spacing Header 

J4 2 Pin .1" Spacing Header 

J5 3 Pin .1" Spacing Header 

J6 6 Pin .1" Spacing Header 

J7 DB9 Right Angle Through Hole Connector 

R1 100K Pot Through Hole" 

R2 1 OK Thiough Hole 

R3 10K Through Hole 

R4 330 Through Hole 

R5 330 Through Hole 

SW1 Momentary Push Button Switch SW403-ND 
SW2 Momentary Push Button Switch SW403-ND 
U1 PIC 18F4620 (40pin DIP) 

U2 74HC14 (14pin DIP) 

U3 AD558 (16pin DIP) 

U4 lO.OOOMliz Crystal Oscillator Full Size 

U5 7405 5V Regulator (T0-220) 

U6 MAX3221 (16pin SSOP) 
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Frequency Converter PIC18F4620 C Code: 


//Jacob Osterberg 
//Spring 2009 
//Frequency Converter 
//fieqMain.c 

//♦♦Remember to set the configuration bits 
//♦♦♦♦Oscillator — HS-PLL enabled freq=4xFoscl 

//****PQRTB A/D Enable — PORTB<4:0> configured as digital I/O on RESET 
//♦♦♦♦£ 0W voltage Program -- Disabled 

#include <pl8t4620.h> 

#include "frequency.h" 

void lowisr(void): 
void highisr(void); 

#pragma code high_vector=0x08 
void interrupt_at_high_vector(void) 

I 

asm GOTO highisr endasm 

} 

#pragma code /* return to the default code section */ 

#pragma interrupt highisr 
void high_isr (void) 

{ 

if(INTCONbits.INTOIF) { 

FREQ_ISR(): 

} 

} 

#pragma code low_vectoi-=0xl8 
void intemipt_at_low_vector(void) 

{ 

asm GOTO low isr endasm 

} 

#pragma code /* return to the default code section */ 

#pragma iutemiptlow low_isr 
void low isr (void) 

{ 

if(INT C ONbits. TMROIF) { 

TimerOISRO; 

/ 

} 

void main() { 
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Function_InitO; 

INTCONbits.GIE = 1: //global enable 

TIMER O RUN = 1: //start the timer 

while(l){ 

FreqToDAC(); 

} 

} 


//Frequency.h 

#ifiidef_FREQ 

#defme_FREQ 

#incl«de <delays.h> 

#include <timers.h> 

#iuclude <math.h> 

#define HIGH 1 
#defme LOW 0 
#define INPUT 1 
#define OUTPUT 0 
#defme TRUE 1 
#define FALSE 0 
#defme ON 1 
#defme OFF 0 

#define TIMER O RUN TOCONbits.TMROON 
#defme TIMER O FLAG INTCONbits.TMROIF 

#defme TMRO_CONST 2500000. // 1/2500000 = seconds per timer increment 

//Constants 

#defme MAX FREQ 713 //Freq Value (Hz) resulting in maximum DAC output (<2500Hz) 

#define MIN FREQ 509 //Freq Value (Hz) resulting in minimum DAC output (>40Hz) 

#defme FREQ SCALE 60 //Set largest with MAX_FREQ*FREQ_SCALE < 65535 

#defme DATA OUTPUT COUNT 14 //# of values averaged before output - (between 1-30) 

#defme DAC_MAX VAL 255 //# of DAC Values 

//****************************************************************** 

#defme MAX_VAL (unsigned) MAX_FREQ*FREQ_SCALE 
#defme MIN_VAL (unsigned) MIN_FREQ*FREQ_SCALE 
#define FREQ STEP (unsigned) ((MAX_FREQ*FREQ_SCALE)- 
(MIN_FREQ*FREQ_SCALE))/DAC_MAX_VAL 

void Function_lnit0; 
void FreqToDAC'(): 
void FREQ_ISR(): 

void Timer0_ISR(): 
void Timer0_lnit(): 
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void WriteDAC(char); 
void Delay 100: 

#endif 


//Frequency.c 

#include "firequency.h" 

unsigned int FREQ_L = 0: 
unsigned int FREQ_H = 0; 
unsigned int FREQ = 0: 
unsigned int ave FREQ = 0: 
unsigned long totalFREQ = 0: 

char DataFlag = FALSE: 
char OverflowFlag = FALSE: 
char DataCount = 0: 

//Initialization Function for interrupts, inputs, outputs and timerO 
void Function_Init0{ 

//Interrupts 

RCONbits.IPEN = ON: //enable interrupt priority levels 
INTCONbits.PEIE = ON: //peripheral interrupt enable 

//Setup INTO 

TRISBbits.TRISBO = INPUT: //RB0 INTO is an input 
INTCON2bits.INTEDGO = HIGH: //trigger interrupt on rising edge 
INTCONbits.INTOIF = FALSE: //clear the INTO flag 
INTCONbits.INTOIE = ON: //enable INTO 

//setup DAC OUTPUT 

TRISD = OUTPUT: //set D port as output 

PORTD = 0: //initialize as 0 

TRISCbits.TRISC3 = OUTPUT: //set C3 as output for DAC latch 
PORTCbits.RC3 = 0: //initialize as 0 


//Initialize Timer Zero 

Timer0_lnit0; 

} 

//Initialization of timerO 
void TimerO_IuitO{ 

//timer 0 register setup 

//16 bit internal clock timer with 2x scalar 

//initially turned off 

TOCON = ObOOOOOOO 1: 

INTCONbits.TMROIE = ON: //enable timer 0 interrupt 

INTCONbits.TMROIF = FALSE; //make sine the time 0 overflow flag not set 
INTCON2bits.TMROIP = LOW: //low priority interrupt 

//set timer to 0 
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TMROH = OxOO: 

TMROL = 0x00: 

} 

//Sets the overflow flag if the timer overflows 
void TimerO_ISR(){ 

INTCONbits.INTOIE = OFF: //disable INTO - don’t let it harware interrupt while 

//resetting the timer 

//notify system of overflow 
Overflow_Flag = TRUF: 

//clear the interrupt flag 
INTCONbits.TMROIF = FAFSE: 

//stop the timer 
TIMER_0_RUN = OFF: 

//enable INTO 

INTCONbits.INTOIE = ON: 

} 

//Takes an 8 bit value and writes it to the DAC 
void Write_DAC(char data){ 

PORTD = data; //write data 

Delayl0(): 

PORTCbits.RC3 = 0: //open latch for write 

Delayl0(): 

PORTCbits.RCS = 1: //latch in data 

} 

//Rims when interrupt tiggers - gets the value of the timer, averages over DATAOUTPUTCOUNT 
//values, resets the timer, checks for overflow, and signals for output after 
//DATA OUTPUT COUNT values 
void FREQISR0 { 

unsigned int Time_L: 

if(Overflow_Flag = FALSE) { 

TIMER_0_RUN = OFF; //stop the timer 

//read the timer (have to read low byte first) 

Time_F = TMR0F; 

FREQ = TMROH: 

FREQ «= 8: 

FREQ = TimeL: 

//Reset Timer 
TMROH = 0x00: 

TMROL = 0x00: 

//run the timer 
TIMER 0 RUN = ON: 


126 



//These conditional statements only signal valid data for output 
//(by setting the Data_Flag) when the number of measured frequency values 
//is greater than DATA OUTPUT COUNT <-set hi frequency.il 
if(Data_Couut < DAT AOUTPUTCOUNT) { 
totalFREQ += FREQ: 

Data Coimt++; 


} 

else{ 


} 


aveFREQ = total_FREQ/DATA 
totalFREQ = 0; 

Data_Couut = 0; 

DataFlag = TRUE: 


OUTPUTCOUNT: 

//reset the frequency sum 
//reset the data count 

//set flag to output data 


//reset the interrupt 

INTC ONbits .INT OIF = FALSE: 

} 

else{ 

//reset timer 
TMROH = 0x00: 

TMR0L = 0x00: 

//TIMER OVERFLOW OUTPUT 
Write_DAC(0): 

DataFlag = FALSE: 

//reset the freq intermpt 
INTCONbits.INTOIF = FALSE: 
//start the timer again 
TIMERORUN = ON: 

//clear the overflow 
OverflowFlag = FALSE: 


} 

//Watches for new data and outputs to the DAC when data is available 
void FreqToDAC(){ 

long value = 0: 
unsigned data = 0: 


//if there is new data 
if(Data_Flag) { 

//find the scaled frequency value 

value = TMR0_CONST*FREQ_SCALE/ave_FREQ: 

//see if the freqeucy is less than the allowed minimum 
//if so output 0 to the DAC 
if(value <= (MIN_VAL)){ 

WriteDAC(O); 

} 

//if the frequency is greater than the maximum value output 
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//set the output to the maximum DAC output value 
else if(value >= (MAX_VAL)){ 

WriteD AC (DAC _M AX_ V AL): 

} 

//if the frequency is within the correct range output it to the DAC 
else{ 


data = (value - MIN_VAL)/((tmsigned) FREQ STEP): 
if(data < 255 && data > 0) 

Write_DAC((char) data): 

} 

DataFlag = FALSE: //reset the new data flag 

/ 


} 


//Delays for 10 cycles 
void Delayl0(){ 

Delay 1TCY() 
Delay 1TCY() 
DelaylTCY() 
DelaylTCY() 
DelaylTCY() 
Delay 1TCY() 
DelaylTCY() 
Delay 1TCY() 
DelaylTCY() 
DelaylTCYQ 
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